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Abstract. We present the results of a multiband survey for brown dwarfs in the Chamaeleon I dark cloud with 
the Wide Field Imager (WFI) camera at the ESO/MPG 2.2-m telescope on La Silla (Chile). The survey has 
revealed a substantial population of brown dwarfs in this southern star forming region. Candidates were selected 
from R, I and Ha imaging observations. We also observed in two medium-band filters, M855 and M915, for the 
purpose of spectral type determination. The former filter covers a wavelength range containing spectral features 
characteristic of M-dwarfs, while the latter lies in a relatively featureless wavelength region for these late-type 
objects. A correlation was found between spectral type and (M855-M915) colour index for mid- to late M-type 
objects and early L-type dwarfs. With this method, we identify most of our object candidates as being of spectral 
type M5 or later. Our results show that there is no strong drop in the number of objects for the latest spectral 
types, hence brown dwarfs may be as abundant as low-mass stars in this region. Also, both kind of objects have a 
similar spatial distribution. We derive an index a — 0.6 ±0.1 of the mass function in this region of dispersed star 
formation, in good agreement with the values obtained in other star forming regions and young clusters. Some 
of the brown dwarfs have strong Hq emission, suggesting mass accretion. For objects with published infrared 
photometry, we find that strong Ha emission is related to a mid-infrared excess, indicative of the existence of a 
circumstellar disk. 

Key words, stars: low-mass, brown dwarfs - stars: pre-main sequence - stars: formation - stars: luminosity function, 
mass function - stars: circumstellar matter 
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1. Introduction vary from one region to another. For instance, some Orion 

clusters are rich in substellar objects (Lucas & Roche OUUUl 

Our understanding of substellar objects has considerably aimB, while in Taurus only a few brown dwarfs 

increased in the last years. Recent surveys have revealed ^^^^^ (Briceiio et al.IMI Martin et allMTO. 

important populations of brown dwarfs in young open ^^^^^^ environmental conditions may play an important 

clusters and starforming regions (e^ouvier et al.|MI ^ole in the formation of brown dwarfs. Stellar winds from 

B^r et al. Wilking et al. M Lucas & Roche ^ ^^^^^^ ^^^^-^^ ^^^^ ^^^1^ evaporate the envelopes 



Barrado y Navascues et al. IMIl Zapatero Osorio accreting cores before they reach a mass above the hy- 

et al.EDEB Lamm et al. in prep.). Still, little is known ^^^^^^ ^^^^^^ ^.^^^ j^^^^^^^ possibility is that brown 

about the properties of very young substellar objects (j-^^rfs are stellar embryos whose growth stopped after 

^ they were dynamically ejected from small stellar systems 

An open question is the way brown dwarfs form. (Reipurth & ClarkeEDUIl Bate et al. Brown dwarfs 

Although it is commonly accepted that they originate, like ^^^^y also form like planets, in accretion disks around stars 

low-mass stars, from the gravitational collapse of molecu- g_ pj^^g^t et al. POOO ) . In this case, most of them would 

lar cloud cores, their formation process is still poorly un- j^^^g ^een ejected from their parental systems to account 

derstood. The number of observed brown dwarfs seems to so-called "brown dwarf desert" (e.g. Sterzik & 

~ ; ~ ~~ ~~ -n T ' Durisen 1999). Another way to explain the lack of brown 

bend offprint requests to: a. Lopez Marti i r i i i 

* r, ' J V . . u 4- J 4- 4-1, T? o 4-1, dwarf companions m wide orbits around stars is that they 

Based on observations collected at the European Southern ^ ■' 

Observatory La Silla Chile have migrated inwards due to their high mass (compared 
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to planets). Eventually, they could remain in closer orbits, 
or continue their migration until they merge with the cen- 
tral star (Armitage & Bonnell lWI^ . 

The Chamaeleon I cloud is a favourable region for the 
study of the formation of substellar objects, because of its 
proximity (~160 pc, Whittet et al. I1997|l and high galac- 
tic latitude {b « —16°), which decreases the effects of con- 
tamination by other late-type stars in the Galaxy. It is the 
largest of the three dark clouds in the Chamaeleon com- 
plex (Gregorio Hetem et al. 11989(1 . and apparently also the 
oldest. The cloud contains more than 150 known young 
stars (e.g. Gauvin & Strom Feigelson et al. [T?M1 

Appenzeller & Mundt I1989|l with a mean age of about 
3 Myr, most of them clustered in two cloud cores contain- 
ing two intermediate mass stars, HD 97048 and HD 97300. 
Chamaeleon I has also been surveyed for brown dwarfs: 
Neuhauser & Comeron H1998|l reported the discovery of 
the very first X-ray emitting brown dwarf, ChaHal, in 
Chamaeleon L In subsequent observations, another twelve 
bona-fide brown dwarfs and brown dwarf candidates (with 
masses between about 0.08 and O.OSMq) were found in 
a small region around HD 97048 (Comeron, Rieke & 
Neuhauser 1999; Comeron, Neuhauser & Kaas 2000; here- 
after EEMni and 101101 respectively). 

This paper presents a new large and deep survey for 
brown dwarfs in the Chamaeleon I cloud. In Sect.|21we de- 
scribe the reduction and photometry of our data. Our can- 
didate selection criteria follow in Sect. 01 A method for the 
spectral type identification of M and early L-dwarfs from 
direct imaging is presented. We discuss our object clas- 
sification, and compare our identifications with available 
near- and mid-infrared data. We then proceed to study the 
possible binarity of the low-mass objects in Chamaeleon I 
(Sect.El, their spatial distribution (Sect.^, their accre- 
tion indicators (Sect.lHJl and their mass function (Sect.O. 
The implications of our results for the proposed brown 
dwarf formation models are discussed in Sect.|Hl We sum- 
marize our conclusions in Sect. El 

2. Observations and Data Processing 

2.1. WFI Observations 

The Wide Field Imager (WFI) is a focal reducer-type 
camera mounted at the Cassegrain focus of the 2.2-m 
ESO/MPG telescope at La Silla (Baade et al. [TM^ . It 
consists of a mosaic of eight 2 x 4 K CCD detectors with 
narrow interchip gaps (the filling factor is 95.9%). The de- 
tector has a field of view of about 34' x 33' and a pixel 
size of 0'.'24. 

In Chamaeleon I, our survey consists of four WFI 
fields, covering an area of about 1.2 (see Fig.QJ. Two of 
these fields contain the intermediate mass stars HD 97048 
and HD 97300, respectively, around which most of the 
young low-mass stellar population is clustered. The field 
previously surveyed for brown dwarfs jCkN99l KJNi^OOll 
covered a region around HD 97048 that represents about 
15% of the total area of our ChaI-4 field. The other two 



fields lie in regions where no intermediate mass stars are 
present: ChaI-5 to the south of ChaI-4, and ChaI-6 in the 
region between ChaI-4 and ChaI-7. We did not dithered 
between the fields, so some objects might be missed be- 
cause they fall in an interchip gap. 

The observations were carried out on three different 
nights in May/June 1999. We observed in the broad-band 
filters R and I, for the purpose of candidate selection, and 
in a Ha filter, to test the youth of our objects. Because 
we wanted to attempt a photometric spectral classifica- 
tion (see Sect. l3.2ll . we also observed in two medium-band 
filters, 855/20 and 915/28^ (hereafter M855 and M915, 
respectively), centred inside and outside TiO and VO ab- 
sorption bands. These bands get deeper with later M- 
subspectral type, but disappear in K- and L-objects. 

For each field and filter, we took three exposures with 
different exposure times, in order to prevent the bright- 
est T Tauri stars (TTS) from being saturated. For the 
Ha filter only two different exposure times were used. For 
the field ChaT7 this observing programme could not be 
carried out completely due to technical problems, and in 
this case we only have images with short and intermediate 
exposure times in the R and I filters. The average seeing 
varied from - 1" for the field ChaL7 to - 1'.'7 for ChaI-5. 
The log of our observations is shown in Tabled 

Bias exposures and skyflats were taken in each fil- 
ter. For the broad-band photometry calibration, several 
Landolt fields (Landolt |1992j) were observed in all nights. 
We also observed a list of stars with known late spectral 
type (between K7 and L3) to use them as standards for 
the spectral type calibration (see Sect. I3.2|l . 

2.2. Data Reduction and Photometry 

The image processing (data reduction and photometry) 
was run with the MSCRED and CCDRED packages in 
the IRAF environment.^ Due to the large field of view, a 
standard reduction of WFI images (bias subtraction and 
flat field division) is not enough to obtain good photome- 
try results: First, because the illumination of the sky flat 
fields and the science images is not equivalent, there re- 
mains an illumination gradient in the background of our 
object exposures after the flatfielding. To correct for this 
effect, the image has to be divided through this gradient 
or illumination mask. Second, due to the long exposure 
times, a fringing pattern appears. This effect is especially 
strong in the I band. To correct for it, one has to subtract 
an individually scaled fringe mask from each science im- 
age. The fringes cannot be completely removed, however, 

^ The transmission of these filters was remeasured on 
December 16th 1999, and they were re-named 856/14 and 
914/27, respectively (see the WFI web site for details). 
However, we will keep the old nomenclature throughout this 
paper. 

^ IRAF is distributed by the National Optical Astronomy 
Observatory (NOAO), which isoperated by the Association of 
Universities for Research in Astronomy, Inc. under contract to 
the National Science Foundation. 
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Fig. 1. DSS image of the Chamaeleon I cloud showing 
the location of our four WFI fields, ChaI-4 through 7. The 
open circles indicate the positions of the objects studied by 
ICRN99l and lCNK00l the filled circles and triangles, of our 
candidates with and without Ha emission, respectively. 

Table 1. Log of WFI observations in Chamaeleon I 



Date 


Field 


Filter 


Exp. Times (s) 


28 May 1999 


ChaL4 


Rc/162 


600/60/5 


28 May 1999 


ChaL4 


Ic/lwp 


600/30/5 


28 May 1999 


ChaL4 


Hq/7 


600/15 


28 May 1999 


ChaL4 


855/20 


600/300/16 


28 May 1999 


ChaL4 


915/28 


600/100/8 


28 May 1999 


ChaL6 


Rc/162 


600/60/5 


28 May 1999 


ChaL6 


Ic/lwp 


600/30/5 


28 May 1999 


ChaL6 


Ha/7 


600/15 


28 May 1999 


ChaL6 


855/20 


600/300/16 


28 May 1999 


ChaL6 


915/28 


600/100/ 8 


30 May 1999 


ChaL7 


Rc/162 


60/5 


30 May 1999 


ChaL7 


Ic/lwp 


30/5 


30 May 1999 


ChaL7 


Hq/7 


600/15 


30 May 1999 


ChaL7 


855/20 


600/300/16 


30 May 1999 


ChaL7 


915/28 


600/100/8 


3 Jun 1999 


ChaL5 


Rc/162 


600/60/5 


3 Jun 1999 


ChaL5 


Ic/lwp 


600/30/5 


3 Jun 1999 


ChaL5 


Hq/7 


600/15 


3 Jun 1999 


ChaL5 


855/20 


600/300/16 


3 Jun 1999 


ChaL5 


915/28 


600/100/8 



because their intensity varies from one chip to the other 
and even between different areas of the same chip. On the 
other hand, in the R and Ha images, as well as in most of 
the images with short exposure times, the fringing effect 
is negligible, and therefore we only corrected for the il- 
lumination gradient in those cases. Both the illumination 



and the fringe masks were created by combination of the 
science exposures. 

For the astrometry, we followed the standard method 
proposed in Francisco Valdes' Guide to the NO AO Mosaic 
Data Handling Software.^ Stars from the Guide Star 
Catalogue were selected as reference for the astrometric 
calibration (Morrison et al. 121)1)11 and references therein). 
The astrometric error is in general not larger than 1". 

Object search was then performed on each individual 
chip of the long exposures using SExtractor (Bertin & 
Arnouts [l99()|l . After several tests, we concluded that this 
software worked better than DAOFIND, leading to fewer 
spurious detections. 

For the photometry, the DAOPHOT package was used 
(Stetson TWI. Since the DAOPHOT tasks only work on 
individual CCDs, a special IRAF procedure was written 
for the WFI images. We first performed aperture photome- 
try with the PHOT task, and then PSF-photometry using 
the PSF and ALLSTAR tasks, regardless of the crowd- 
edness of the field. In this way we intended to measure 
the faintest components of multiple systems. This did not 
prevent some very faint objects to be merged with their 
much brighter neighbours, however, although we allowed 
only for a very small merging radius. This problem affects 
to objects up to separations of about 4" depending on the 
relative brightness of the objects. Given that the fraction 
of faint objects lost due to this problem is small (less than 
3%) in such an uncrowded region, we decided not to fur- 
ther consider these objects in our study. Only in a few 
interesting cases (see Sect. EJ, aperture photometry was 
performed on the subtracted images. 

On most of the short exposures only aperture photom- 
etry was performed, because of the difficulty of fitting a 
suitable PSF. This was also the case of the standard stars 
(exposure times of 5-10 s). 

2.3. Photometric Calibration 

The two broad-band filters were calibrated with Landolt 
1)1992(1 standard stars. Absolute magnitudes were com- 
puted using the relations: 

R = r + Bi+B2-XR + B3-{r~i) (1) 

I = t + Ai+A2-Xj + A3-{r-i) (2) 

where r and i are the instrumental magnitudes, and 
Xj denote the measured airmasses, and Ak and Bk are the 
coefficients obtained from a linear fit using the standard 
instrumental and absolute magnitudes. 

In the first two nights (May 28th and May 30th), five 
Landolt fields containing more than 60 standard stars were 
observed. Their airmasses range from 1.15 to 2.90. The fits 
for these nights have a standard deviation of about 5%. 
In the fourth night (June 3rd), only two Landolt fields 
were taken, although one of them was observed twice at 

^ This guide is available from: 
,http: //iraf .noao . edu/ scripts/iraf help?mscguide] 
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Table 2. Calibration coefficients for the WFI I and R 
filters. 



Night 


Al 


A2 


Ai 


Bi 


Sa 




May 28th 


-1.38 


-0.13 


0.23 


-0.21 


-0.24 


0.00 


May 30th 


-1.47 


-0.03 


0.25 


-0.58 


0.00 


0.00 


June 3rd 


-1.19 


-0.19 


0.32 


-0.18 


-0.23 


0.00 



different airmass. The number of stars used for the fit (the 
stars in the repeated field were used twice) was 29, with 
airmasses between 1.14 and 1.95. The standard deviation 
is about 6%. 

The resulting coefficients are somewhat different for 
each night (see Table In particular the colour term for 
the I band, A3, is relatively high. This is not surprising, 
because the shape of the transmission curve for this filter 
differs notably from a standard Cousins I filter, especially 
at shorter wavelengths. On the other hand, the R filter 
shows a good match with a Cousins R filter, and in this 
case the colour term is so small that it can be set to zero 
(see the WFI User Manual for details). 

The main error source in this calibration is the fact 
that it comes from a global fit for all the CCD chips. 
Recently, Alcala et al. H2UU2|I showed that, for the WFI, 
chip-to-chip photometric variations can be as high as 3% 
in the broad-band filters and 5% in medium-band filters, 
because of different zero-points for the individual CCDs. 
With no doubt, a better result could be obtained by fitting 
the stars in each frame separately. This was not possible 
in our case, however, because some CCD chips contained 
few or no standard stars (especially those in the corners). 
Future observations should take several exposures dither- 
ing the fields, so that enough stars are observed on each 
CCD. 

We also note that the Landolt stars used for the photo- 
metric calibration span a range in (R~I) colours of about 
-0.1 to 1.5, while our objects of interest have colours in the 
range 0.5 to 4 (see Sect. I3.ljl . Hence, our colours may not 
be exactly in the Cousins system, especially for the red- 
dest objects with the largest colour terms. This is a general 
problem that affects all studies of very red objects, as long 
as no set of very red standard stars is available. For the 
WFI colours, comparison with the predictions of the Lyon 
models (Baraffe et al. HMHl Chabrier et al. IMTHjl shows 
deviations of up to 0.5 mag for very red objects objects 
(R-I=2.5), while there is good agreement for R-I=1.5 (see 
also Eisloffel et al., in prep.). 

Apart from this, our photometry shows a systematic 
error of unknown origin: Using the objects present in the 
images with different exposure times to estimate the error 
in our photometry, we found a systematic offset between 
them. In general, the resulting luminosities are brighter 
by a constant amount (typically between 0. and 0.06 mag, 
though as high as ^ 0.1 mag in the worst case) for longer 
exposure times. This effect is reproduced with various 



photometry packages and with different parameters. It is 
present already in the raw data (i.e., previous to the re- 
duction) , and does not show any obvious correlation with 
the observed field, the filter used or the measured extinc- 
tion. Observing conditions are not likely to have caused 
it, because the images were taken one after another, and 
the effect is always in the same direction for all images in 
all photometric nights. Extensive tests by the ESO staff 
have verified that the effect is not caused by a WFI shut- 
ter malfunction leading to systematically short or vari- 
able integration times. It seems also improbable that it is 
caused by detector saturation, because the objects appear 
brighter, not fainter, with longer exposure times. 

Therefore, to homogenize our photometry, we have de- 
termined the offsets between the three different integra- 
tion times, and have shifted the short and long integration 
times to the system of the medium integration. That way, 
we should have minimized this systematic error, since the 
medium integration times for our objects are not very dif- 
ferent from those of the standard stars. We did not choose 
the shortest exposures as reference because they contained 
only a few objects which could be used to estimate the ex- 
posure time offsets. 

Since our observed fields in Chamaeleon I are slightly 
overlapping, we checked for double detections of objects. 
We found that the photometry results in the overlapping 
region between two fields are consistent within the esti- 
mated errors. Comparison with DENIS I photometry for 
previously known objects in our field ChaI-4 (Neuhauser & 
Comer6n ll999ll gives differences of about 0.2 mag. Hence, 
our results are marginally consistent with the DENIS mea- 
surements, whose errors for faint objects can get as high 
as 0.2 mag for I~18 mag (Epchtcin et al. 11997(1 . A larger 
difference is found when comparing our RI photometry 
with the results of l(:;RN99l and l(:;NK00l Our results yield 
systematically fainter magnitudes by about 0.4 mag in 
average in both passbands. This is not surprising, given 
that the photometry of lCRN99l and lCNK00l was calibrated 
with few secondary standards in the observed field, which 
are probably variable, and without a colour correction. 
Because our error seems to be approximately equal for 
both the R and I photometry, our R— I colours should be 
very little affected by the exposure time effect. 

Because no standard stars are available for the 
medium-band filters, no absolute calibration was possible 
for the Ha, M855 and M915 photometry, and only instru- 
mental magnitudes will be considered for the data in these 
filters. We applied a standard correction for atmospheric 
extinction, using the extinction coefficients of the R and 
I band for the Ha and the two M filters, respectively. For 
a better understanding of the emission properties, the Ha 
magnitudes were then shifted so that the main locus of the 
objects in each surveyed field corresponds to Ha— R=0 in 
a (Ha, Ha~R) colour- magnitude diagram (see Sect. |3| be- 
low). 

To estimate the completeness limits of our survey, 
we constructed the number distribution of the found ob- 
jects as a function of their measured magnitudes for our 
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Table 3. Photometry and spectral types of the well known objects in Chamaeleon l^^ 



5 



— KT 


a (2UUUJ 


(2U0Uj 


R 


I 


Ha 


M855 


M915 


— 7. G 

SpT 


SpT 


Classification 




11 08 54.7 


-77 32 12.3 


17.50 


15.12 


17.23 


18.92 


18.22 


M5.5 


M5.5 


star 




CHJvH 21 


11 07 11.8 


-77 46 37.8 


15.28 


13.46 


15.04 


17.24 


16.86 


K7 


<M4 


star 




CHXK 22 


11 07 13.7 


-77 43 50.3 


17.21 


14.96 


16.94 


18.65 


18. 14 


"^1 TV T 1^ 

K7+JV12 


M4 


star 




CHXR 74 


11 06 57.8 


-77 42 11.2 


15.92 


13.91 


15.59 


17.70 


17.20 


M4.5 


M4 


star 




HM 19 


11 08 ib.y 


-77 44 37.8 


15.12 


13.31 


14.88 


17.18 


16.76 


M3.5 


<M4 


star 




Sz 23 


11 07 58.5 


-77 42 42.0 


16.55 


14.64 


15.59 


18.43 


18.03 


M2.5 


<M4 


star 




ChaHa 1 


11 07 17.2 


-77 35 53.5 


19.45 


16.62 


18.95 


20.54 


19.40 


M7.5 


M9 


BD 




ChaHcK 2 


11 07 42.9 


-77 34 00.2 


17.96 


15.42 


17. 15 


19.19 


18.44 


M6.5 


M6 


trans. 


obj. 


V^llcLxlCt O 




77 "ifi 1^7 ft 


17 80 


15 26 


17 61 


19 15 


IS 30 


M7 


M6 5 


trans. 


obj. 


ChaHa 4 


11 08 19.4 


-77 39 17.8 


17.08 


14.74 


16.85 


18.63 


17.89 


M6 


M6 


trans. 


obj. 


ChaHa 5 


11 08 24.6 


-77 41 48.1 


17.56 


15.02 


17.36 


18.89 


18.06 


M6 


M6.5 


trans. 


obj. 


ChaHa 8 


11 07 46.5 


-77 40 09.6 


18.37 


15.78 


18.22 


19.67 


18.79 


M6.5 


M7 


BD 




ChaHa 9 


11 07 19.2 


-77 32 52.4 


20.6; 


17.65 


19.95 


21.30 


20.46 


M6 


M6.5 


trans. 


obj. 


ChaHa 10 


11 08 24.5 


-77 39 30.7 


19.86 


17.25 


19.74 


21.09 


20.22 


M7.5 


M7 


BD 




ChaHa 11 


11 08 29.7 


-77 39 20.4 


20.4; 


17.72 


19.89 


21.59 


20.65 


M8 


M7 


BD 




ChaHa 12 


11 06 38.4 


-77 43 09.7 


18.78 


16.09 


18.45 


19.96 


18.96 


M7 


M8 


BD 




ChaHa 13 


11 08 17.5 


-77 44 12.5 


16.95 


14.52 


16.64 


18.41 


17.59 


M5 


M6.5 


trans. 


obj. 



^ Hq, M855 and M915 magnitudes arc instrumental magnitudes (sec text). 

^ Photometric errors: blank: .05; semicolon: 0.1 mag; colon: 0.2 mag 

Spectral types from lCNKOOl fb rown dwarfs and brown dwarf candidates) and ICRN99I (stars'!. 

Spectral types from this work. The estimated errors are 1 subclass in the range M6-M9 and 2 subclasses in the range M4-M6 (see text). 



four WFI fields. These distributions peak for values of 
R~20 mag and I~19 mag, except for the field ChaI-4, 
where these values are shifted <1 mag towards fainter ob- 
jects for both filters. There is a sharp drop in the num- 
ber of objects for R~23 mag (22 mag for ChaI-7, due to 
the shorter exposure times) and I~21 mag, respectively. 
Inspection of the corresponding images shows that objects 
fainter than these values are most likely noise fluctua- 
tions. Hence, we consider general detection limits of about 
23 mag and 21 mag for the R and I filters, respectively, 
and completeness limits of R~20 mag and I~19 mag, al- 
though the survey might be complete for fainter objects 
in the area covered by the field ChaI-4. 

3. Survey Results 

3.1. Brown Dwarf Candidates 

At an early stage of their evolution, brown dwarf candi- 
dates in star forming regions are easily identified in (I, 
R-I) colour-magnitude diagrams. Fig. |21 shows our dia- 
gram for Chamaeleon I. In the already surveyed region 
around the star HD 97048, we recover all the previously 
known objects from IURN991 and lUNEOOl ( except for two 
of them that happened to fall in the interchip gaps): They 
lie in the region I<14 and are clearly separated from back- 
ground objects towards redder (R~I) colours. Our optical 
photometry for these objects is summarized in Table |21 
These previously known objects have an estimated age 
of about 1-3 Myr, and define an empirical isochrone in 
the colour-magnitude diagram, which can be used to se- 
lect candidates in a much larger field. Such an isochrone 
contains information about the distance and reddening of 
the objects in the dark cloud. Moreover, by using an em- 
pirical rather than a theoretical isochrone, we avoid the 
uncertainties still present in the evolutionary models for 
very young ages (see Baraffe et al. 120021 for a critical dis- 
cussion) . 




12 3 4 
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Fig. 2. (I, R-I) colour-magnitude diagram for our sur- 
veyed field in Chamaeleon I. The previously identified low- 
mass stars and brown dwarfs from CRN99 and CN KOOl 
(open symbols) define an empirical isochrone towards red- 
der colours, where we select our new candidates (solid 
symbols). The crossed lines indicate the average errors. 



Using such an isochrone, we find 113 candidates mem- 
bers of Chamaeleon I in the same region of the diagram 
and its prolongation towards fainter sources. The solid line 
in Fig. 121 is a fit on the brown dwarfs and brown dwarf 
candidates from ICNKOOL shifted to allow for the spread 
observed in the previously known objects and for the esti- 
mated photometric errors. Our selected candidates cover 
the magnitude range down to I ~ 20 mag. The objects 
tend to be clustered around the cloud cores (see Fig. 
The region most devoid of candidate objects lies towards 
the south of the surveyed area, where only one object is 
found in the zone of interest of the diagram. 

Given the proximity of Chamaeleon I (^160 pc), we 
do not expect a significant contamination from foreground 
objects among our candidates. Moreover, due to its high 
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<, 1-6 - 



1.2 - 




-0.8 -0.4 

Hd-R 

Fig. 3. Logarithm of the Ha equivalent width versus the 
(Ha-R) colour index for the objects in Chamaeleon I 
with measured equivalent width. The data are taken from 
[CRN99 for the low-mass stars, and from CNKOO. for the 
brown dwarfs and brown dwarf candidates. 

galactic latitude (6 ~ —16) and to extinction through the 
dark cloud, we should not be able either to see many back- 
ground objects. To estimate the amount of contaminating 
objects in our field, we made use of the simulations of stel- 
lar population synthesis in the Galaxy by Robin & Creze 
P986|l '*. Within a solid angle of 1 square degree in the 
direction of the Chamaeleon Complex, these simulations 
yield a number of 5 to 7 objects in the region of the (I, 
R-I) colour-magnitude diagram around a 1 Myr isochrone 
from the evolutionary models of Chabrier et al. H2()()()|l . 
Since we find more than hundred objects in the same re- 
gion of the diagram, we are confident that about 95% of 
our candidates are true members of Chamaeleon L 

We note that our survey also includes a number of 
detections in the Tband that have no counterpart in R. 
Most of them lie near the lower magnitude limit of our 
survey with values of I = 20-21. These objects may be 
very red, but because of their non-detection in R nothing 
can be said about their nature. 

According to the theoretical models (Baraffe et al. 
119981 Chabrier et al. I2(JU(J|I . if our candidates are indeed 
members of the star forming region, they should reach 
down to very low masses, extending down to the deu- 
terium burning limit for the faintest objects in our sample. 
One way to confirm cloud membership is the detection of 
Ha emission, which is usually considered as an indicator 
of youth. All the objects from lCR,N99l and ICNKOOl indeed 
show the Ha line in emission in their spectra. 

The amount of emission is usually quantified by the 
measurement of the Ha equivalent width EW(Ha). Fig.O 
shows a plot of logEW(Ha) versus our measured (Ha-R) 
colour index for the previously known brown dwarf and 
brown dwarf candidates. The values of EW(Ha) are taken 
from ICNKOOl With the notable exception of ChaHal, 



16 



20 




-3 -2 -1 



Ha-R 

Fig. 4. (Ha, Ha-R) colour-magnitude diagram for our 
surveyed field in Chamaeleon I. Symbols as in Fig. 13 
Objects with Ha emission are placed in the left part of 
the diagram. 

that will be later discussed in Sect. 16.11 a tendency of 
smaller equivalent width towards higher colour index can 
be clearly seen. Of course, one has to be careful with such 
estimations, because the equivalent width measurements 
of ICNKO Ol were performed at an epoch different from our 
optical observations, and these objects are known to have 
variable Ha emission. However, still the (Ha-R) colour 
can be used as a rough estimate of the Ha emission. 

In Fig. ^ we present a (Ha, Ha-R) colour-magnitude 
diagram of the objects in our field. The objects from 
l(^RN99l and ICNKOOl a,re found in the left part of the di- 
agram, having Ha-R colours less than -0.1. Most of our 
bright candidates are also located in this region of the di- 
agram, and thus have Ha in emission. Unfortunately, this 
cannot be clearly stated for most of the faintest objects, 
due to the large photometric errors in the lower part of our 
diagram. However, some of them have such extreme mea- 
sured (Ha-R) colours that they are likely Ha emitters. 
Moreover, many of them also show an infrared excess (see 
Sect 13. 3|) . another indicator that they are probably young. 
Thus, all the candidates with Ha-R<-0.1 in Fig.Qlare con- 
sidered probable Chamaeleon 1 members. In this way, we 
retain 69 objects. Of course, some of the candidates with- 
out apparent Ha emission will also belong to the cloud, 
but this cannot be proved without spectra. 

We later added three other objects to this Hst: 
Chal 608, in spite of being sHghtly too blue for our initial 
selection, is the probable optical counterpart of the X-ray 
source CHXR 30 (Feigelson et al- HM^ . Chal 404 has no 
clear Ha emission, but it could be the optical counterpart 
of an ISOCAM source (see Sect. I3.3|) . Chal 616 has a re- 
ported infrared excess (see Sect. l3.3|l . The coordinates and 
cross-identifications of all 72 objects are listed in Tabled 



* The simulations can be generated on-line at this URL: 
http : //www. obs-besancon. fr/www/modele/modele_aiig. html 
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Table 4. Candidate low-mass members of Chamaeleon I' 



Name 


a 


(2000) 


<5 (2000) 


R 


I 


Hq 


M855 


M915 


SpT'^ 


Classification 


Other Names'^ 


Chal 401 


11 


11 


59.3 


-77 


30 


38.2 


22.9: 


19.9; 


21.3: 


23.2: 


22.9: 


LI. 5 


BD caiid. 




Chal 403 


11 


10 


54.1 


-77 


25 


01.2 


21.0; 


18.05 


19.6; 


21.66 


20.83 


M6.5 


trans, obj. 


DENIS-P JlllO. 9-7725 


Chal 404° 


11 


10 


37.0 


-77 


22 


13.8 


19.70 


16.62 


19.7; 


20.23 


19.36 


M6.5 


BD cand. 


ISO 250 *' 


ChaT 405 




10 


35.4 


-77 


22 


05.9 


18 33 


15 60 


18 16 


19 22 


18 53 


M5 5 




Tqn 9'^n ^ 


Chal 406 


11 


10 


42.0 


-77 


20 


48.8 


20.06 


17.27 


18.25 


21.06 


20.06 


M8 (M2) 


star- 


ISO 252 


Chal 408 


11 


10 


51.4 


-77 


18 


04.0 


16.65 


14.50 


16.41 


18.19 


17.70 


M4 


star 




Chal 410 


11 


08 


19.0 


-77 


30 


41.4 


19.42 


16.90 


18.99 


20.74 


19.85 


M7 (M5.5) 


BD cand. 


ISO 1 3H 


Chal 411 


11 


08 


22.9 


-77 


30 


28.4 


18.02 


15.65 


17.01 


19.42 


18.75 


M5 




KG 45 ISO 143 


PViaT 41 9 


11 


09 


53.9 


-77 


28 


37.2 


21 4' 


18 26 


20 33 


21 98 


20 87 


MQ ('M4 5^1 


star 


PK" "HD mo 990 


PViaT 41 7 


11 


09 


48.0 


-77 


26 


29.8 


18 27 


16 00 


17 07 


19 61 


19 12 


M4 


star 


RaiiH 4"^ Kd 1 01 
mo 907 

low iiVJ 1 


Chal 418 


11 


09 


43.1 


-77 


25 


58.4 


18.90 


15.97 


18.17 


19.61 


18.78 


M6.5 


trans obj 


P7-1 KP q5 ISO 200 


ChaT 41 Q 


11 


10 


02.5 


-77 


25 


45.9 


21 0" 


17 74 


20 4' 


21 32 


20 44 


M7 


BD 




ChaT 4-90 


11 


08 


13.0 


-77 


19 


13.4 


22 1' 


19 03 


20 7' 


22 5' 


21 72 


MQ Q '^ I'M! 1 




KP "^7 mo 1 "^1 
-ivvj oi , low lOl, 

DENIS-P J1108. 2-7719 


Chal 422 


11 


03 


42.5 


-77 


26 


52.8 


18.57 


15.96 


18.19 


19.75 


18.87 


M7 


BD 


ISO 28 


Chal 423 


11 


02 


42.5 


-77 


24 


25.2 


17.71 


15.34 


17.41 


19.20 


18.47 


M6 


trans obj 




PViaT 494 


11 


02 


54.5 


-77 


22 


57.0 


22 5' 


19 3" 


21 2" 


23 3' 


22 1' 


M8 5 


BD 




Chal 425 


11 


03 


48.3 


-77 


19 


57.4 


16.67 


14.22 


16.39 


18.12 


17.20 


M7.5 


BD 


Hn2, ISO 32, DENIS 19 


Chal 427 


11 


03 


11.7 


-77 


36 


36.7 


22.0: 


18.72 


21.0: 


22.4- 


21.48 


M9 


BD 


DENIS-P J1103. 2-7736 


Chal 428 


11 


02 


47.4 


-77 


38 


08.5 


22.6: 


19.0' 


21.2: 


22.4- 


21.74 


LO (Ml) 




DENIS-P J1102. 8-7738 


Chal 429 


11 


07 


24.8 


-77 


43 


49.4 


19.77 


16.93 


19.64 


20.67 


19.77 


M7 


BD 


KG 10 




11 


06 


15.9 


-77 


37 


50.5 


19 07 


16 46 


18 92 


20 02 


19 50 


M4 




ISO 76, 

DENIS-P J1106. 3-7737 




11 


06 


29.3 


-77 


37 


33.6 


18 26 


15 91 


18 09 


19 56 


19 05 


M4 




PHX'T^7'^ ISO 7S 

Wll-^vlV ( O , low (O 


Chal 434 


11 


09 


45.6 


-77 


40 


34.0 


18.02 


15.37 


17.75 


19.22 


18.37 


M7 


BD 


KP 96 ISO 201 




11 


11 


22.8 


-77 


45 


43.6 


21 8" 


18 75 


20 0' 


22 6" 


2 1 50 


M8.5 


BD 






11 


12 


04.0 


-77 


26 


01.5 


18 51 


16 29 


17 85 


20 08 


19 46 


M5 




iqn 989 

lOv^ 


Chal 441 


11 


09 


01.7 


-77 


20 


53.1 


21.0: 


18.42 


20.8: 


22.1- 


21.63 


LI 


BD 


ISO 171 


Chal 448 


11 


07 


21.0 


-77 


29 


41.0 


15.10 


13.22 


14.86 


17.11 


16.59 


M4 




ISO 99 


Chal 449 


11 


05 


43.6 


-77 


26 


52.3 


15.85 


13.65 


15.61 


17.54 


16.84 


M6 


trans obj 


CHXR15 ISO 65 


Chal 450 


11 


06 


44.1 


-77 


26 


35.0 


14.87 


13.11 


14.69 


16.94 


16.56 


<M4 




T22 UX Cha, ISO 80 
DENIS 38 


PhaT 4^9 


11 


06 


59.7 


-77 


18 


54.3 


14.74 


13 02 


13 90 


16 72 


16 49 


^M4 I'M! 




T9'^ TTV Pha 


Chal 453 


11 


02 


33.3 


-77 


29 


13.6 


14.72 


13.11 


14.49 


17.03 


16.69 


<M4 (M3) 


star 


CHXR71, Hnl, 




























ISO 4 DENIS 16 


Chal 456 


11 


02 


55.7 


-77 


21 


51.5 


15.17 


13.40 


14.56 


17.35 


16.77 


M5 (MO. 5) 


star 


T12. SzlO DENIS 17 


Chal 457 


11 


04 


43.1 


-77 


41 


57.7 


15.82 


13.96 


15.49 


17.82 


17.36 


M4 


star 


B18, DENIS 26 


Chal 601 


11 


12 


31.6 


-76 


53 


35.2 


19.35 


16.76 


19.12 


20.75 


19.72 


M8 


BD 




Chal 602 


11 


09 


29.7 


-76 


59 


18.6 


18.04 


15.79 


17.85 


19.67 


19.01 


M5 


star 




Chal 603 


11 


04 


57.6 


-77 


15 


57.2 


16.65 


14.74 


16.43 


18.51 


18.05 


M4 


star 


T16, Szl3, ISO 55, 
DENIS 29 


Chal 607 


11 


08 


12.1 


-77 


18 


54.0 


22.8: 


19.9; 


20.4; 


23.2: 


22.6; 


LO (MO. 5) 


star 


KG36, ISO 130, 
DENIS-P J1108. 2-7718 


Chal 608" 


11 


08 


00.6 


-77 


17 


31.2 


17.79 


15.51 


17.82 


19.26 


18.72 


M4 


star 


CHXR 30, KG 30, 



DENIS 58 



3.2. Spectral Types 

3.2.1. A Photometric Classification 

To confirm the membership of our candidates to the 
Chamaeleon I cloud, we should determine their spectral 
type. Since there are many objects spread over a large re- 
gion, and therefore spectroscopy is not easy, we attempt 
to identify late-type objects exclusively from photometric 
observations. 

M-dwarfs are usually classified by means of their TiO 
and VO absorption features, that get deeper with later 
M-subspectral type (Kirkpatrick et al. 119911 Martin et al. 
I1999|l . but disappear in K- and L-type objects. While TiO 
bands begin to saturate for objects of spectral type later 
than M5, VO features become very prominent for late M- 
type objects. The WFI has a medium-band filter, M855, 
covering some important TiO and VO absorption features 
that get deeper with later spectral type in M-objects (see 
Fig. 5 in Kirkpatrick et al. I1991|l . Another medium-band 
filter, M915, lies in a wavelength range which is relatively 
featureless in late-type objects. For the calibration, sev- 
eral objects with known spectral type from K7 to L3 were 



observed, and a correlation was established between the 
spectral type and the (M855-M915) colour index. 

We selected our spectral calibrators mostly from the 
Gliese-Jahreiss catalogue (H. Jahreiss, private commu- 
nication). Suitable objects must be observable from the 
southern hemisphere, not form multiple systems, and not 
be variable. Unfortunately, the catalogue does not con- 
tain many stars of the latest spectral types (M8-L3), and 
very few fulfilled the required conditions. For this rea- 
son, we also selected some objects from the LHS catalogue 
(Luytcn 1979) and from the surveys DENIS (Delfosse et 
al. 1999 Martin et al. 1999 ) and 2MASS (Kirkpatrick et 
al. I21)UI)|I . although in these cases we could not be sure 
that the objects were single-star systems. 

A first set of calibrators was observed in two nights 
(on the 30th and 31st of May 1999) of the observing run 
in which our Chamaeleon I fields were taken. Once the 
feasibility of our method was demonstrated, a second and 
larger set was observed on the 23rd December 2001 to 
improve our initial calibration. In total, we observed 22 
objects in the M855 and M915 filters, but three of them 
were rejected because they seemed double, or because of 



Chal 731 
Chal 735 



Chal 737 
Chal 742 
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Table |3] (continued) 



Name 


a 


(2000) 


<5 (2000) 


R 


I 


Hq 


M855 


M915 


SpT" 


Classification 


Other Namcs^* 


Chal 609 


11 


07 


57, 


8 


-77 


17 


26.9 


19.78 


17.27 


19.3; 


20.91 


20.32 


M5 


star 


Baud 38, KG 26, 
































DENIS 55 


Chal 610 


11 


08 


00. 


.5 


-77 


15 


32.4 


18.09 


15.49 


17.89 


19.34 


18.47 


M7 


BD 


KG 29, ISO 121 


Chal 613 


11 


08 


27. 


,1 


-77 


15 


55.8 


19.38 


16.73 


18.84 


20.60 


19.62 


MB 


BD 


ISO 147 


Chal 615 


11 


09 


05. 


7 


-77 


09 


59.0 


16.02 


14.09 


15.76 


17.99 


17.45 


M4 


star 


Hn7, ISO 174, 
































DENIS 71 


Chal 616° 


11 


09 


58. 


.8 


-76 


59 


15.5 


14.78 


13.22 


14.76 


17.11 


16.78 


<M4 


star 


KG 113, ISO 229 


Chal 618 


11 


06 


50. 


,4 


-76 


59 


29.1 


16.97 


15.03 


16.86 


20.15 


20.00 


<M4 


star 




Chal 619 


11 


04 


37. 


2 


-77 


02 


14.5 


13.26 


11.89 


13.12 


15.77 


15.53 


<M4 


star 




Chal 622 


11 


05 


15. 


2 


-77 


11 


29.7 


15.10 


13.26 


14.82 


17.09 


16.65 


M4 


star 


Hn4, DENIS 30 


Chal 625 


11 


10 


29. 


,1 


-77 


17 


00.4 


16.18 


14.04 


15.89 


17.92 


17.21 


M6 


trans, obj. 


Hnl2, DENIS 94 


Chal 629 


11 


11 


23. 


.3 


-77 


05 


54.7 


15.87 


13.95 


15.64 


17.80 


17.36 


<M4 


star 




Chal 708 


11 


09 


13. 


.9 


-76 


28 


40.1 


15.89 


13.99 


15.61 


18.03 


17.43 


M5 


star 


Hn 8, KG 75, 
































DENIS 73 


Chal 709 


11 


08 


51, 


,0 


-76 


25 


13.8 


16.81 


14.79 


16.26 


18.87 


18.20 


M5.5 


star 


T37, Sz28, KG 54, 
































ISO 157, DENIS 68 


Chal 710 


11 


10 


22. 


,3 


-76 


25 


14.2 


19.09 


16.54 


18.62 


20.66 


19.50 


M9 


BD 




Chal 711 


11 


07 


11. 


,9 


-76 


25 


50.2 


18.37 


16.15 


18.06 


20.01 


19.32 


M5.5 


star 




Chal 712 


11 


06 


32. 


,8 


-76 


25 


21.0 


19.42 


16.96 


19.1; 


20.93 


19.92 


M8 


BD 




Chal 717 


11 


08 


02. 


,5 


-76 


40 


34.3 


17.92 


15.54 


17.68 


19.58 


18.68 


M7 


BD 




Chal 721 


11 


09 


51. 


,6 


-76 


45 


45.4 


16.06 


14.30 


15.92 


18.09 


17.78 


<M4 


star 


ISO 214 


Chal 726 


11 


09 


52. 


,2 


-76 


39 


13.0 


19.59 


16.96 


18.45 


20.87 


19.92 


M7.5 


BD 


GK 29, KG 106, 



11 09 22.8 
11 08 55.1 



11 08 51.9 
11 12 03.4 



-76 34 32.2 
-76 32 40.1 



-76 32 50.3 
-76 37 03.7 



19.94 
18.03 



20.3; 
16.28 



17.55 
15.76 



17.61 
14.13 



18.88 
17.12 



20.0; 
15.90 



21.11 
19.63 



21.60 
18.18 



20.74 
18.95 



20.62 
17.43 



<M4 
M5.5 



MB 
M6 



star 
star 

BD 

trans, obj. 



ISO 217 

Cl-6, NIR 10, ISO 1B9 
Cl-25, NIR 19, KG 59, 
ISO 165 

Hnl6, CHXR 84, 
DENIS 105 



Chal 


743 


11 


11 


45.4 


-76 


36 


50.1 


21.4: 


IB, 


.60 


21.0: 


22.8; 


21.47 


MB 


BD 




Chal 


747 


11 


13 


24.6 


-76 


29 


23.2 


15.29 


13, 


.69 


14.94 


17.67 


17.32 


<M4 


star 


HnlB, DENIS 116 


Chal 


748 


11 


13 


29.8 


-76 


29 


01.7 


15.23 


13, 


.51 


14.94 


17.48 


17.06 


<M4 


star 


Hnl9, DENIS 118 


Chal 


749 


11 


12 


59.4 


-76 


16 


53.6 


13.57 


11 


.93 


13.23 


15.98 


15.63 


<M4 


star 




Chal 


751 


11 


06 


42.0 


-76 


35 


49.0 


15.28 


13, 


.63 


14.45 


17.61 


17.13 


M4 


star 




Chal 


752 


11 


08 


40.8 


-76 


36 


07. B 


14.15 


12, 


.68 


13.77 


16.65 


16.38 


<M4 (M2) 


star 


CHX 13a, KG 52, 
ISO 153, DENIS 65 


Chal 


753 


11 


09 


53.5 


-76 


34 


25.5 


15.10 


13, 


.52 


13.87 


17.31 


17.09 


<M4 


star 




Chal 


754 


11 


12 


48.8 


-76 


47 


06.7 


15.46 


13, 


.79 


15.15 


17.78 


17.34 


<M4 


star 




Chal 


755 


11 


10 


56.2 


-76 


45 


32.7 


16.28 


13, 


.93 


15.67 


17.94 


17.02 


M7 


BD 




Chal 


757 


11 


12 


10.0 


-76 


34 


36.8 


15.00 


13, 


.12 


14.51 


17.18 


16.60 


M4 


star 




Chal 


759 


11 


13 


33.7 


-76 


35 


37. B 


15.37 


13, 


.59 


15.08 


17.57 


17.11 


M4 


star 




Chal 


760 


11 


10 


11.5 


-76 


35 


29.3 


15.56 


13, 


.89 


15.29 


17.74 


17.47 


<M4 (Ml. 5) 


star 


NIR 45, KG 121, 



ISO 237, 

DENIS-P JlllO. 2-7635 



Hq, MB55 and M915 magnitudes are instrumental magnitudes (see text). 
^ Photometric errors: blank: 0.05; semicolon: 0.1 mag; colon: 0.2 mag 

° The estimated errors in the spectral types arc 1 subclass in the range M6-M9 and 2 subclasses in the range M4-M6 (sec text). Spectral types 
from the literature are indicated in parenthesis. 

References: Sz#: S chwar tz <l<)77t : Hn#: Harti gan JltiOiSI : T#: Schwartz iliiiiH; CHX R#: Feigelson et al. y^jictt: C#: Prusti et al. M§11: 

NIR#: Oasa et al. J10i>{)l : Baud#: Baud et al. Jljiij41 ; CHX# : Feigelson & Kriss llOijOl B#: Lawson e t al. il^gSl : ISOi: Persi et al. ii^OoSl . P. 

Persi, priv. communication; DEN IS#: C ambresy et al. I199BI Table 1); DENIS-P#: Cambresy et al. 119981 'I'able 2); GK#: Gomez & Kenyon 

12223; KG#: Kenyon & Gomez 120011 
^ This object docs not have Ha emission according to its (Ha— R) colour. 
^ Ambiguous counterpart (see text). 

^ The (R-I) colour of this object is slightly too blue according to our selection criteria. 



technical problems during the observations. Data reduc- 
tion and aperture photometry were performed in the way 
explained in Sect 12.21 Tabled summarizes our photometry 
results and other relevant data of our final calibrator list. 
We recall that no absolute photometric system is avail- 
able for the M855 and M915 passbands. We could not use 
these objects to establish an absolute magnitude scale in 
these filters because for most of them the distance is still 
unknown. 

For mid- to late M-type objects (M4-M9), a correla- 
tion between colour and spectral type is clearly found (see 
Fig.O: Objects with later spectral type also have a higher 
(M855-M915) colour index. This is not seen for early 
M-type objects, because the features considered are not 
prominent in their spectra. The observed scatter probably 
is due mainly to differences in age (see also Sect. 13.2. 2|l : 



metallicity differences may also play a role. The best fit 
obtained in this range is: 

SpT = (6.2 ± 1.2) • (M855 - M915) + (2.8 ± 0.7) (3) 

where SpT is an integer (4-9) denoting the M subspectral 
class. The quoted errors are the standard deviations of the 
corresponding coefficients. 

For early L-type objects (L0-L2), an anticorrelation is 
observed, as expected from the progressive weakening of 
the TiO and VO features. However, this anticorrelation 
seems to break down for objects of spectral type later 
than L2. This is probably due to the appearance of the al- 
kali absorption features, a characteristic of the L spectral 
sequence. One Cs I line lies towards the centre of the wave- 
length range covered by our M855 filter, and it starts to be 
important for L3 and later-type objects (see Kirkpatrick 
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IN <xllil_. 




2000) 


S (2000) 


M855 


M915 


SpT 




09 


28 


51 


-09 


16 


00 


13 


12 


13 


19 


K7 


GJ 7 59 


06 


51 


42 


-43 


53 


12 


12 


96 


13 


02 


MO 


GJ 494 


13 


00 


47 


12 


22 


33 


11 


83 


11 


69 


MO 5 


GJ 7 58 


06 


39 


50 


-61 


28 


48 


12 


71 


12 


77 


Ml 


VjJ ZZUOD 


08 


16 


08 


01 


18 


12 


12 


24 


12 


26 


M2 


P T 4 fif^iQ 
vjj 4: Doy 


07 


38 


41 


-21 


13 


30 


13 


59 


13 


56 


M3 


GJ 402 


10 


50 


52 


06 


48 


29 


12 


97 


12 


64 


M3 




05 


36 


00 


-07 


39 


00 


14 


40 


14 


20 


M4 


P T 9 1 (T^i^ 


03 


09 


00 


10 


01 


24 


15 


99 


15 


69 


M5 


GJ 4 195 


02 


02 


16 


10 


20 


18 


16 


19 


15 


73 


M6 


LHS 3003 


14 


56 


38 


-28 


09 


50 


16 


57 


15 


72 


M7 


GJ 5 14 


02 


52 


26 


00 


56 


18 


19 


81 


19 


27 


M8 


LHS 2065 


08 


53 


36 


-03 


29 


32 


18 


48 


17 


40 


M9 


GJ 4 619 


07 


07 


53 


-49 


00 


48 


20 


22 


19 


42 


M9 


DENIS-P J0909. 9-0658 


09 


09 


57 


-06 


58 


06 


20 


67 


20 


19 


LO 


2MASSW J0832045-012835 


08 


32 


05 


-01 


28 


35 


20 


98 


20 


57 


LI. 5 


2MASS J1342236+175156 


13 


42 


36 


17 


51 


56 


19 


81 


19 


80 


L2 


2MASSW J0928397-16031 


09 


28 


40 


-16 


03 


12 


21 


30 


21 


34 


L2 


DENIS-P J1058. 7-1548 


10 


58 


47 


-15 


48 


00 


21 


18 


20 


71 


L2.5 


DENIS-P J1047. 5-1815 


10 


47 


31 


-18 


15 


58 


21 


40 


20 


60 


L2.5 



^ MS55 and M915 magnitudes are instrumental magnitudes. 

^ Photometric errors: blanlt: 0.05; semicolon: 0.1 mag; colon: 0.2 mag 

References: GJ#: H. Ja hreiss (private communication); LHS#: Luyten 119791 : DENIS: Delfosso ct al. Cl999j . Martin et al. 119991 : 2MASS: 

Kirkpatrick et al. 120001 . 



M855-M915 



Fig. 5. Spectral type versus (M855-M915) colour index 
for our set of spectral calibrators. Numbers from to 9 
correspond to the M subspectral classes, and from 10 on 
to the L spectral sequence. For mid- to late M-type objects 
(M4-M9), a correlation between colour and spectral type 
is clearly seen. The bar at the lower right corner of the 
plot shows the average error in the (M855~M915) colours. 



et al. 119991 their Figs. 4, 6 and 7). Since L3 is the hmit 
of our observed spectral range, we cannot state whether a 
sequence similar to that of late M-objects, but based on 
the alkali absorption features, could be established also for 
late L-objects. In the range M9-L2, we obtain as best fit: 

SpT = (-3.1 ± 0.6) • (M855 - M915) + (12.0 ± 0.4) (4) 

where SpT is now in the range 9-12 (a value of 10 corre- 
sponds to a spectral type LO). 



MB55-Mgi& 



Fig. 6. Spectral type versus (M855-M915) colour for the 
objects in Chamaeleon I with known spectral type. The er- 
rors are taken from '^GNKQQI A correlation between colour 
and spectral type is clearly seen. The long-dashed line is 
the fit obtained with our spectral calibrators. The solid 
line is the fit for this second set of objects. Both lines 
are almost parallel, as indicated by the short-dashed line, 
which corresponds to the fit on the old objects, but shifted 
2 subclasses. The dotted lines indicate the estimated error 
of about 1 subclass. 



3.2.2. Application to Young Objects 

In a next step, we used Eq. Oto derive spectral types for 
the known objects in Chamaeleon I. 

The results are shown in Fig. El The young objects 
appear to have, according to our fit, later spectral types 
than the ones derived from spectroscopic observations. 
However, the trend to increasing colour with later spec- 
tral type is clearly reproduced. The solid line shows a fit 



10 



Lopez Marti et al.: The brown dwarf population in the Chamaeleon I cloud 



performed on this set of young objects. This Une is al- 
most parallel to the one obtained from the fit to the (old) 
field objects, just being shifted to earfier spectral types by 
about 2 subclasses. 

We interpret this shift as a consequence of the lower 
value of log 5 for these young objects. Work by several 
authors (e.g. Martin et al. I199f)l Luhman et al. I1997|l 
shows that young M-type objects (characterized by low- 
gravity atmospheres) display stronger TiO and VO ab- 
sorption than field dwarfs of the same spectral class 
(high-gravity atmospheres) . Stronger TiO and VO absorp- 
tion implies fainter M855 magnitudes and redder (M855- 
M915) colours. The age diff'erence should also be respon- 
sible, to some extent, for the scatter of our calibrators. 

The effect of age (i.e. gravity) in our medium-band 
photometry can be tested with the help of theoretical 
models. The left panel of Fig. |7| shows different isochrones 
in a (Te//, M855-M915) diagram from the non-grey mod- 
els of Baraffe et al. (1998; I. Baraffe, private communica- 
tion). The represented ages range from 1 to 600 Myr. All 
curves tend to converge for effective temperatures higher 
than 4000 K, corresponding to late K-type objects. For 
lower temperatures, the younger objects tend to exhibit 
redder (M855-M915) colours than the older ones of same 
spectral type. The biggest differences are observed in the 
temperature range between 4000 and 2800 K, roughly 
corresponding to spectral types from K7 to M6. Here, a 
1 Myr old object will have the same colour as a 600 Myr 
old M-dwarf that is about 800 K cooler. Even compared 
to a 3 Myr old object a remarkable temperature differ- 
ence is predicted 200 K). However, for cooler objects 
(Te// < 2800 K), the colour difference is no longer sig- 
nificant at very young ages (less than 10 Myr), and the 
difference between young and old objects is notably de- 
creased (from 800 to about 150 K). On the other hand, 
the isochrones for intermediate ages (10-70 Myr) appear 
more or less parallel through the whole represented range. 
The oldest isochrone (600 Myr) only diverges from them 
at very low temperatures. 

The relation between spectral type and (M855-M915) 
colour index can be better understood with the help of 
the right panel of Fig.d Here we have adopted the Teff 
scale versus spectral type for pre-main sequence objects 
from Luhman l|1999|l : 

Te// = 3850 - 141.0 •S'pT (5) 

It is immediately seen that the behaviour of the curves is 
not exactly linear: For blue (M855-M915) colours (earlier 
spectral types) the difference between old and young ob- 
jects can be as high as 5 subclasses, but this difference is 
reduced towards redder colours (later spectral types). If we 
shift the oldest isochrone 2 subclasses, however, the differ- 
ence is not greater than 1 subclass for spectral types later 
than M4. The value of this shift is the approximated offset 
that we find between our fit to (old) field M-dwarfs and the 
(young) objects in Chamaeleon I. We therefore conclude 
that the observed offset between our initial spectral type 



calibration and the reported spectral types for the objects 
in Chamaeleon I is mainly due to the difference in age be- 
tween them and our calibrators, which implies a different 
value of logg. To account for this effect, we shifted our 
derived scale 2 subclasses towards earlier spectral types 
(see Fig. O. In most of the cases, the so derived spectral 
types coincide with the ones in ICRN99I and .CNKOOl (ob- 
tained from object spectra) with an error of one subclass 
at most (see TableOJ. For early M-type objects, the error 
is certainly larger, due to the saturation of the (M855- 
M915) colour. We conservatively consider an error of two 
subclasses in the range M4-M6. 

We cannot test our (M855-M915) colour index for L- 
dwarfs in the same way, because no L-type objects are 
yet known in Chamaeleon I. The model curves of Fig. [3 
do not go deep enough into the L-type temperature range 
to test the observed behaviour. Therefore, we can only 
tentatively assign LO to L2 spectral types to the objects 
found at the bottom left part of a (M855, M855-M915) 
colour-magnitude diagram. We note, however, that highly 
extincted M-type objects would also fall in the same re- 
gion of the diagram, making this classification somewhat 
uncertain. Hence, it is difficult, at the present stage, to dis- 
tinguish between a L-type object and a rather extincted 
M-type object with this method (see also the discussion in 
the following section). A more detailed analysis is required 
before a reliable determination of the L subspectral class 
can be achieved with this photometric technique. 

3.2.3. Spectral Types of the Brown Dwarf Candidates 

Fig. IHl shows the (M915, M855-M915) colour-magnitude 
diagrams for our Chamaeleon I field. Our derived scale 
for the spectral type identification is also shown. Clearly, 
many of our candidates have mid to late spectral type M. 
About one fifth of them have spectral type M7 or later. 
Given that the stellar/substellar boundary at the age of 
Chamaeleon I is set at a spectral type of about M6 by 
the theoretical models, we classify these objects as brown 
dwarfs. There is a smaller group of seven or eight objects 
with spectral types between M5.5 and M7, in the transi- 
tion region from low-mass stars to brown dwarfs. Around 
50 objects have spectral types M5 or earlier, so they are 
very low-mass stars. For the faintest objects, a turning 
back towards lower values of the (M855-M915) colour is 
seen in Fig.|Slat spectral type M9, which may be indicative 
of the beginning of the L-dwarf sequence. If their belong- 
ing to the cloud can be confirmed, these objects might be 
low-mass brown dwarfs. 

The photometry results and the derived spectral type 
for all our candidates are summarized in Table 01 The 
values obtained from Eqs.|31and01have been rounded to 
the closest value in steps of 0.5 subclasses. Although 23 
objects from our survey list had already been proposed 
as cloud members in previous optical and/or X-ray stud- 
ies of Chamaeleon I (see references in Table 0J, many of 
them still lacked a spectral type determination. We have 
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MB55-Hgi5 M855-Mgi5 



Fig. 7. Left panel: Isochrones of effective temperature versus theoretical (M855-M915) colour index from the Baraffe 
et al. H1998|l models. The curves of equal mass are also plotted. For the same value of the colour index, older objects 
have a lower temperature, thus a later spectral type. Hence, by using a fit to old objects, the temperatures and masses 
of young very low-mass objects are overestimated. Right panel: Spectral type versus theoretical (M855~M915) colour 
index for the 1 Myr and the 600 Myr isochrones. At low values of the colour index, older objects are several subclasses 
later than very young objects. This difference decreases to at most 1 subclass for objects later than M4 if the oldest 
isochrone is shifted by 2 subclasses. 



checked our photometrically derived spectral types with 
those given by Lawson et al. and Gauvin & Strom 

l|1992|l for the objects in common with these authors (7). 
In general there is good agreement between their and our 
spectral types, except for the star Sz 10 (Chal 456), to 
which Gauvin & Strom assign a spectral type of MO. 5, 
while our photometric classification yields a spectral type 
of M5. We do not see, however, any other object near 
Chal 456 that could be the one observed by Gauvin & 
Strom. DENIS photometry is also available for this ob- 
ject (Cambresy et al. IDMl see Sect. E31 below). The I 
magnitude measured by these authors is consistent with 
our own result within the estimated errors, while their 
JKs photometry agrees well with that given by Gauvin 
& Strom (1992). We are thus quite confident that our 
identification of Chal 456 with Sz 10 is correct. However, 
the measured I magnitude (13.40 mag) appears too bright 
for an M5 object. The discrepancy could be explained if 
Sz 10 were an unresolved binary. In that case, Gauvin & 
Strom would have observed the spectrum of the primary. 
We would wrongly assign a later spectral type to the un- 
resolved pair, because our spectral type calibration is only 
based on the (M855-M915) colour index (which should be 
redder for a binary object), and not on magnitudes (which 
should be brighter). As a matter of fact, in TableOlwe see 
that CHXR22, a known binary formed by a K7 and a M2 
object that appear unresolved in our images, has an esti- 
mated spectral type of M4 according to our classification. 
Thus an hypothetical secondary of Sz 10 would probably 
have a spectral type earlier than M5, definitively being a 
star. 

More recently, Gomez & Persi l)2002|l presented in- 
frared spectroscopic observations of several mid-infrared 
sources from a previous ISOCAM survey in Chamaeleon I 



(Persi et al. I2000|l . Six of their objects are also included 
in our list of probable cloud members (some of them in 
common with other authors) . They give an error of 2 sub- 
classes for their spectral type classification. While for four 
of the objects the agreement between their and our spec- 
tral types is rather good, for the other two we find again a 
large discrepancy: ISO 252 (our object Chal 406) is clas- 
sified as an M2 object according to these authors, but 
we derive a spectral type of M8 for it. Its I magnitude 
according to our and the DENIS measurement (around 
17.2 mag) is slightly too faint compared to that of other 
objects of similar spectral type in the cloud. On the other 
hand, the object ISO 250 (Chal 405), which lies not too 
far from ISO 252, has a spectral type of M6 according to 
Gomez & Persi and of M5.5 according to our own classi- 
fication, in spite of being brighter than ISO 252 in the I 
band. The discrepancy between both spectral type clas- 
sifications for ISO 252 might be due to a locally higher 
extinction towards this object. A similar situation is seen 
in the case of ISO 220 (Chal 412), which we classify as M9 
and Gomez & Persi as M4.5. The I magnitude for this ob- 
ject (greater than 18 mag according to both DENIS and 
our photometry) is very faint compared to other objects 
of spectral type M4-M5 in the cloud (typically between 
16 and 17 mag). The (R-I) colours of both Chal 406 and 
Chal 412 are too red for early M-type objects, which might 
indeed hint to the presence of moderate extinction. 

We note, however, that the discrepancy between lumi- 
nosity and spectral type might also be due to the pres- 
ence of circumstellar dust, as proposed by Fernandez & 
Comeron (|2001|1 for the object LS-RCrA 1 in Corona 
Australis. This object appears too faint compared to other 
objects of similar spectral type and age. More such sublu- 



12 



Lopez Marti et al.: The brown dwarf population in the Chamaeleon I cloud 



■ ■ m • 

° "(^ l.o 



1 

°* o. 



1 ^ >* ""^ 



MB55-M915 



Fig. 8. (M915, M855-M915) colour-magnitude diagram 
for our surveyed field in Chamaeleon I. Symbols as in 
Fig. 121 A scale for the identification of the spectral type is 
also indicated (see text). 



minous objects have been identified in other star forming 
regions. 

The effect of extinction is more clearly seen when com- 
paring our spectral types with those given by Gomez & 
Mardones 



These authors performed near-infrared 
spectroscopy of 46 candidate young stellar objects from 
several photometric surveys in Chamaeleon I. We have 
six objects in common with this work.^ For three of them 
(Chal 403, Chal 432 and Chal 760), the spectral classifica- 
tion given by Gomez & Mardones agrees with ours within 
their and our estimated errors. This is not the case, how- 
ever, for the other three objects (Chal 420, Chal 428 and 
Chal 607), whose spectral types are remarkably warmer 
according to the spectroscopic data. While we assign them 
spectral types between M9 and LO with our photometric 
method, Gomez &: Mardones classify them as early M-type 
objects (M0.5-M1). These authors obtain values between 
2 and 5 mag for the extinction Aj in the J band towards 
these objects, which are among the faintest ones in our 
survey. 

These cases illustrate the limitations of our photomet- 
ric spectral type determination: We will estimate too late 
spectral types for some objects that lie in regions of locally 
high extinction, which could then be wrongly classified as 
brown dwarfs. The uncertainties are especially large for 
the faintest candidates in our sample, that could be iden- 
tified as L-type objects in spite of being highly extincted 
early M-type stars. 

On the other hand, if a particular object is a bi- 
nary which appears unresolved in our WFI images, we 
are likely to assign a later spectral type to the combined 
pair. Nevertheless, among a total of 36 objects with known 
spectral type, we have only six clear discrepancies, mean- 
ing that in about 80% of the cases our classification is 



^ Gomez & Mardones II2UU31 actually observed seven of our 
objects, but the spectrum of one of them, Chal 427, was too 
noisy to allow a spectral type classification. 



reliable. This reliability is even higher (around 90%) for 
the brightest objects (I<19). On these objects, a linear 
fit of our predicted spectral types versus the ones from 
the literature gives a RMS error of about one subclass, in 
agreement with our expectations. 

Out of 72 objects, we have classified 44 objects as low- 
mass stars (with spectral type earlier than M6), and 19 ob- 
jects as brown dwarfs (with spectral type later than M6.5). 
We also find 6 transition objects with spectral types be- 
tween M6 and M6.5, which will be considered brown dwarf 
candidates in the subsequent discussion. Three objects 
that could be low-mass brown dwarfs of spectral type L 
are also detected. However, in some cases the identifica- 
tion of our faint candidates as cloud members is not cer- 
tain because of the large photometry errors in the lower 
part of the (Ha, Ha-R) colour-magnitude diagram and/or 
the lack of infrared excess (see Sect. l3.3|l . Moreover, as ex- 
plained in this section, their spectral type classification 
may be somewhat uncertain and requires spectroscopic 
confirmation. We will regard these objects as brown dwarf 
candidates, although keeping in mind that they might be 
embedded stars. 

In the subsequent analysis we do not include the 
doubtful objects Chal 608 and Chal 404. In the seven 
cases mentioned above where the photometric spectral 
type does not match with reported spectroscopic obser- 
vations, the classification derived from object spectra was 
preferred. Spectral types from the literature are indicated 
in parenthesis in the corresponding column of Table 01 

3.3. Infrared Detections 

We cross-correlated our Chamaeleon I object list with 
published detections by ISOCAM (Persi et al. l^nnUjl and 
DENIS (Cambresy et al. 11998(1 , as well as the surveys by 
Gomez & Kenyon (20011) and Kenyon & Gomez (2001). In 
this way, we got infrared photometry for about 67% of the 
objects in our sample. Of them, 47 have a near-infrared 
counterpart and 33 a mid-infrared counterpart. 

Persi et al. H2()()()(l provide a list of faint Chamaeleon I 
members with mid-infrared excesses from ISOCAM ob- 
servations. We found 19 objects in common in both lists. 
Another 11 objects from our sample have been detected by 
ISOCAM in one or both bands (P. Persi, private commu- 
nication). The ISOCAM list also provides DENIS IJKs 
photometry for all the detected sources. In addition, we 
recover 23 objects with near-infrared excess from DENIS 
photometry according to Cambresy et al. (19981 . Gomez 
et al. H2001|l report a list of 118 candidate members of the 
Chamaeleon I cloud with infrared excess from JHK ob- 
servations, but we find only two objects in common with 
these authors. Kenyon & Gomez H20()l|l provide JHKsL 
photometry for a sample of young stellar candidates in 
Chal, including many from the previous surveys and some 
new detections. This sample contains 16 of our objects. 
In total, we find 33 objects with a reported near-infrared 
excess. For the ones observed by several authors, the pho- 
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Name 


Jd 


Kso 


J 


H 


K 


L 


1^6.7 (mjy) 


1 14.3 (mJy) 


Chal 403 


13.28 


10.83 
















Onai 4u4 


12.34 


10.80 












20.0 ± 0.0 


5.0 it 1.2 


Chal 405 


12.34 


10.80 












20.0 ± 0.0 


5.0 it 1.2 


Chal 406 


13.75 


12.30 












5.5 ± 0.0 


6.0 it 1.2 


Chal 410 


14.07 


13.02 












15.3 ± 0.0 


17.6 it 0.3 


Chal 411 


12.65 


11.06 


12 


.74 


11.65 


11.10 


10.67 


15.6 ± 1.0 


16.5 it 1.5 


Chal 412^^ 


14.36 


12.38 


14 


.57 


13.33 


12.48 




5.3± 0.0 


6.0 it 1.1 


Chal 417 


12.38 


9.93 


12 


.74 


11.12 


10.30 


8.58 


34.4 ± 1.0 


35.0 it 1.4 


Chal 418 


12.19 


10.42 


12 


.48 


11.12 


10.48 


9.55 


14.9 ± 1.0 


12.3 it 2.3 


Chal 420 


13.10 


8.83 


12 


.96 


10.28 


8.90 


7.74 


13.5 ± 1.0 


36.2 it 1.4 


Chal 422 


13.05 


11.91 












3.2 it 0.0 




Chal 425 


11.26 


9.96 












11.4 ± 0.0 


3.7 it 1.0 


Chal 427 


13.40 


10.01 
















Chal 428 


13.75 


10.50 
















Chal 429 






13 


.62 


12.47 


11.83 


10.72 






Chal 432 


12.72 


10.08 












12.6 ± 1.0 




Chal 433 


12.78 


10.51 












7.0 ± 1.0 




Chal 434 


12.39 


10.93 


11 


.92 


10.92 


10.47 


10.03 


7.4 lb 1.0 




Chal 438 


13.01 


11.56 












7.6 ± 1.0 


8.0 it 1.5 


Chal 441 


14.55 


11.69 












3.4 ± 0.0 




/^l T 4 A O 

Chal 448 


11.30 


10.14 












11.2 ± 0.0 




/^l T 4 An 

Chal 449 


11.29 


10.09 












13.2 ± 1.0 




Chal 450 


10.94 


9.27 












20.7 ± 0.0 




Chal 453 


11.23 


10.16 












15.0 ± 0.0 


7.5 it 1.9 


Chal 456 


11.41 


10.36 
















Chal 457 


11.78 


10.61 
















Chal 603 


12.10 


10.27 












23.1 ± 1.0 


14.4 it 2.1 


Chal 607 


14.54 


10.48 


14 


.46 


11.72 


10.47 


9.92 


14.3 ± 1.0 




Chal 608 


11.74 


9.04 


11 


.71 


9.86 


9.08 


8.45 






Chal 609 


13.08 


9.87 


13 


.19 


11.00 


9.61 


8.29 






Chal 610 


12.31 


11.13 


12, 


.46 


11.57 


11.10 


10.38 


4.8 ± 1.0 




Chal 613 


13.73 


12.34 












4.3 ± 1.0 




Chal 615 


11.79 


10.77 












5.9 ± 0.0 




Chal 616 


10.46 


8.73 


10 


.59 


9.38 


8.94 


8.51 


26.6 it 1.0 


6.0 it 1.2 


Chal 622 


10.93 


9.46 
















Chal 625 


11.60 


10.64 
















Onai 7Uo 


11.64 


10.67 


11 


.81 


11.34 


10.79 


10.28 






Chal 709 


12.45 


11.30 


12 


.38 


11.66 


11.25 


10.53 


7.5 lb 1.0 




Chal 721 


11.10 


9.26 












24.0 lb 1.0 




Chal 726 


13.28 


11.70 


13 


.36 


12.34 


11.65 


10.76 


10.3 lb 2.0 




Chal 731 


13.13 


9.32 












383.1 lb 5.0 


611.3 it 19.3 


Chal 735 


13.03 


11.38 


12 


.92 


11.96 


11.36 


10.90 


12.0 lb 1.0 




Chal 742 


11.69 


10.66 
















Chal 747 


11.78 


10.64 
















Chal 748 


11.53 


10.45 
















Chal 752 


10.45 


9.09 


10 


.56 


9.56 


9.20 


8.92 


26.8 lb 2.0 




Chal 760 


10.77 


8.51 


10 


.87 


9.34 


8.55 


7.69 


186.3 lb 6.0 


238.0 it 8.8 



" JoKsd-- DENIS photometry (Cambrosy et al. 119981 : JHKL: Photometry from Kenyon & Gomez 120011 . 
^ Photometry from Gomez & Kenyon 120011 



tometry results are consistent within the quoted errors. 
The available infrared photometry for our objects, to- 
gether with the cross-identifications and references are 
summarized in Table |H| Surprisingly, only 10 objects ex- 
hibit both near- and mid-infrared excesses. This represents 
less than 25% of the total infrared detections in our sam- 
ple. 

A dubious cross-identification regards the objects 
Chal 404 and Chal 405, which lie relatively close at a 
distance of about 9'.'5. The coordinates given by Persi et 
al. ll^UnUjl for the ISOCAM source ISO 250 match well 
the position of Chal 404, the fainter of the two, but the 1- 
band photometry from DENIS most probably corresponds 
to Chal 405, which is about one magnitude brighter than 
Chal 404. In our "reddest" filter, M915, the brightness dif- 
ference between both sources is reduced to about 0.8 mag. 
It is possible that Chal 404 is brighter in the infrared than 
its neighbour, thus leading to an eventual identification er- 
ror. However, Chal 404 does not show evident Ha emission 



according to its Ha-R colour, so its membership to the 
Chamaeleon I cloud is doubtful according to our selection 
criteria. Given the ISOCAM pointing errors, we cannot 
exclude that Chal 404 is indeed a background object, and 
that Chal 405 is the true counterpart to the mid-infrared 
source. 

The detection of infrared excesses indicates that 
the Chamaeleon I candidate members must be young. 
However, it should be taken into account that criteria for 
identifying objects with infrared excess are not homoge- 
nous among the different authors. In particular, we note 
that the Gomez & Kenyon survey has few objects in com- 
mon with the other two, in spite of their spatial overlap. 
For a detailed discussion of the infrared properties, we re- 
fer to the quoted works. Nonetheless, something has to be 
said about the extinction towards our objects, as it may 
affect their properties, and in particular, according to the 
previous section, their spectral type classification. 
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In Fig. 121 we show the near-infrared colour-colour dia- 
gram for our candidates with JHK photometry (Gomez 
& Kenvon OnUn Kenyon & Gomez "^IMl). The solid curve 
is a 1 Myr isochrone from the Chabrier et al. H2(J(J01|) mod- 
els. The dashed and dotted lines indicate the position of 
this isochrone for extinction values of Ay = 5 mag and 
Ay = 10 mag, respectively, computed with the law of 
Rieke & Lebofsky H1985|l . The colours of most of our can- 
didates are consistent with extinction values between 
and 10 mag. The three exceptions are Chal 420, Chal 607 
and Chal 609, that appear to be significantly reddened. 
For the two former objects, this result was expected after 
the discussion from Sect. 13.231 Note that for the third ob- 
ject with a remarkably discrepant spectral type (Chal 428) 
only DENIS photometry is available. Chal 609, on the 
other hand, is very red for its estimated spectral type 
(M5). 

Persi et al. PTOjl used a {J - K, K-m{6.7)) colour- 
colour diagram to separate the effects from intrinsic in- 
frared excess from those of reddening. 

The apparent magnitudes m(6.7) and m(14.3) corre- 
sponding to the ISOCAM fluxes (in Jy) are computed with 
the following relations: 

m(6.7)^-2.5.1og^:il|^ (6) 

m(14.3).-2.5.1og^:^iiiif^ (7) 

In Fig. ^1 we show the position of our objects in 
the (J — K, K—m{6.7)) diagram. The solid and dashed 
lines indicate the location of the reddening band com- 
puted by Persi et al. (2000). The objects with intrinsic 
(i^T— m(6.7)) excess are found to the right of the reddening 
band. Chal 420 and Chal 607 are clearly separated from 
the rest, having high values of the (J — K) colour index, 
but low or moderate (X— m(6.7)) colours. This diagram 
confirms that these two objects are heavily reddened, and 
that many from our candidates suffer from moderate ex- 
tinction. 

4. Spatial Distribution 

Fig-fflshows the positions of our objects the Chamaeleon I 
region, as well as those of the previously known brown 
dwarfs and very low-mass stars from CRN99, and .CNKQOi 
The distribution of the brown dwarfs and brown dwarf 
candidates resembles that of the low-mass stars. All these 
objects are mostly located in the two cloud cores con- 
taining the two intermediate mass stars HD 97048 and 
HD 97300 and at the boundaries of a third core, appar- 
ently the densest according to C^^O observations (Mizuno 
et al. I1999|l . This core lies between the other two and is 
seen as a dark region devoid of stars in the optical images. 
Far from the three cloud cores, few objects are found. 

To better quantify this observed distribution, wc have 
counted the number of objects at different distance bins 
from the stars HD 97048 and HD 97300. Fig.lTTlshows the 
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Fig. 9. Colour-colour diagram for our candidates with 
near-infrared photometry. The curves indicate the posi- 
tion of a 1 Myr isochrone from Chabrier et al. H2000|l for 
different extinction values. The crossed lines on the upper 
left corner of the plot indicate the average errors according 
to Kenyon & Gomez H2001fl . 
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Fig. 10. (J — K, K—m.{6.7)) colour-colour diagram for 
our candidates detected in the near and mid-infrared. The 
objects with intrinsic {K—m{6.7)) excess are found to the 
right of the assumed reddening band, indicated by the two 
parallel lines. 

histograms derived from these counts.^ We do see that the 
distributions of both stars and brown dwarfs peak near 
the positions of the intermediate mass stars. In general, 
the objects seem to be more clustered near HD 97048 (a 
Herbig Ae/Be star) than near HD 97300 (a likely zero-age 
main sequence A star). A third peak is observed in both 

® Although these counts might be somehow incomplete, due 
to the fact that our studied area is larger in the North-South 
than in the East- West direction, we do not expect a signifi- 
cant contribution from objects beyond the survey boundaries 
because the cloud is mostly covered by our observations (see 
e.g. Mizuno et al. ll999t . Thus, we do not correct for boundary 
effects, which could introduce a new bias in the analysis given 
that such corrections are generally based on the assumption of 
a random distribution beyond the boundaries. 
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Distance [arcmin] 

Fig. 11. Distribution of the low-mass Cha I members ac- 
cording to the distance to the intermediate mass stars 
HD 97048 (upper panel) and HD 97300 (lower panel). The 
blank histograms show the distribution of the very low- 
mass stars; the hashed histograms, of the brown dwarfs 
and faint brown dwarf candidates. The general histograms 
for all these objects are showed in dotted lines. The dis- 
tance to the other intermediate mass star is also indicated. 



cases at the distance corresponding to the boundaries of 
the third cloud core between both stars (at around 30' 
from HD 97048 and 40' of HD 97300). Although some of 
the objects in the bins could actually be very far away 
from this core, the peak is significant because it appears 
when counting from both stars. 

5. Binary Objects 

To look for possible binary pairs among our objects, we 
adopt an upper separation cutoff of 12" (about 1800 AU). 
This value was determined by estimating the average sur- 
face density in our surveyed area. Then the distance was 
computed at which a given object should have a single 
neighbour if there were a uniform distribution on the 
whole field; this distance turned out to be about 12'.'2. 
Previous studies of visual binaries in Chamaeleon I (e.g. 
Reipurth & Zinnecker '19^21 Brandner et al. 19913 ) also 
used an upper limit of 12". Indeed, inspection of our im- 




Fig. 12. Detail of the long I exposure showing the possi- 
ble multiple system formed by the star Chal 405 (upper 
object), the brown dwarf candidate Chal 404 (lower ob- 
ject) and a very faint object towards the north-east of 
Chal 404. The two brightest components lie at a separa- 
tion of about 9'.'5. This system is found at the position of 
the ISOCAM source ISO 250. (North is up and east is left 
in this image.) 



ages shows that the probability of finding a neighbouring 
object is notably increased at a slightly larger separation. 

Most of our objects are found in isolation, or at least 
do not have neighbours down to the distances resolvable 
from our images. In very few cases one of our objects is 
found near a previously known cloud member: Chal 625, 
classified as a transition object (spectral type M6), is iden- 
tified as the brightest component of the visual binary Hnl2 
fSchwartz lTMT]! . The M4 star Chal 408 is located at about 
14" of the reported position for the star Sz 37 (which lies 
beyond our survey boundaries), and is probably not phys- 
ically bound to it. 

As mentioned in the previous section, the pair formed 
by Chal 404 and Chal 405 lie whithin a separation of 
^ 9'.'5. If it were a true cloud member, Chal 404 would be 
a brown dwarf candidate with an estimated spectral type 
of M6.5. Chal 405 has spectral type M5.5 and is most 
probably a very low-mass star. This system could still 
have a third member at about 5" (~750 AU) of Chal 404. 
This object is only visible in our longest I exposure, with 
a measured magnitude of I~19 mag. Considering our de- 
tection limits, it thus seems to have a very red colour, 
R-I>3, which would fit the empirical isochrone for the 
Chamaeleon members shown in Fig. |2] If this faint ob- 
ject indeed belonged to the cloud, it could be a very low- 
mass companion (a very low-mass brown dwarf or a giant 
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planet) of this possible binary system. We show this pos- 
sible multiple system in Fig. 

Another wide pair with a separation of ^10" is formed 
by Chal 608 and Chal 609. The latter is a star of spectral 
type M5 according to our classification. Its location corre- 
sponds to the infrared-detected star Baud 38 (Baud et al. 
lEHlIl- On the other hand, the position of the ROSAT X- 
ray source CHXR 30 (Feigelson et al. ,1993, ) is well matched 
with that of Chal 608. As mentioned in Sect. 13. ll this ob- 
ject had a slightly too blue (R~I) colour to be included 
in our original candidate selection. We have found only 
one previous reference for a double object in a recent 
paper (Carpenter et al. I2U02|I . though already Feigelson 
et al. reported an offset of 7" between the positions of 
the infrared and the X-ray sources. However, in our own 
inspection of the ROSAT data (Lopez Marti, Stelzer & 
Neuhauser 2002 ) , we do find a bright elongated source at 
the position of this binary candidate, meaning that the 
X-ray source could actually be double and include both 
Chal 608 and Chal 609. 

The pair formed by Chal 735 and Chal 737 (spec- 
tral types M5.5 and M8, respectively), if actually bound, 
would be a good candidate for a wide star-brown dwarf 
system. However, these objects lie at a separation of 
~ 14'.'5, so they have a high probability of not being phys- 
ically connected. 

We also checked our and the previously known objects 
for faint or very close companions that could have been 
missed by the automatic object search. We only find a pos- 
sible companion near a brown dwarf candidate, Chal 425 
(spectral type M7.5), only visible in the longest I expo- 
sure as a faint "tail" to this object. It was not possible to 
get the photometry of this possible companion from the 
subtracted image, due to the residuals left by the subtrac- 
tion of the brightest object (which is almost saturated in 
the longest exposure). On the other hand, as explained in 
Sect. 13. 2. Chal 456 is probably an unresolved stellar bi- 
nary. Chal 726 (a brown dwarf candidate of spectral type 
MS. 5) seems to be slightly elongated in the north-south 
direction and could be a very close binary as well.^ 

Although in several cases faint objects are seen near 
our identified cloud members, their (R~I) colours (when 
both I and R photometry are available) seem to indicate 
that they belong to the background. Also some of the pre- 
viously known brown dwarf candidates have faint neigh- 
bours visible in our images. These objects have recently 
been studied by Neuhauser et al. H2002fl . who conclude 
from their infrared colours that they are unlikely to be- 
long to the Chamaeleon I cloud, and thus to be bound to 
the brown dwarf candidates. To confirm if all these objects 
have close companions, spectroscopic observations as well 
as more imaging observations are needed. 



^ Recalling the discussion from Sect. 13.2.31 if this were the 
case, the individual components of this binary system would 
probably have earlier spectral types than estimated from our 
study, and they might not be brown dwarfs. 
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Fig. 13. Measured (Ha-R) colour index versus spectral 
type. Symbols as in Fig. [21 The strongest emitters are 
among the objects with the latest spectral types. 

From this analysis we conclude that most of our ob- 
jects are not members of (wide) binary systems. Hence, bi- 
narity (down to the separations that we can prove) should 
not significantly change the Chamaeleon I mass function. 

6. Accretion Processes 

6.1. Ha Emission 

The detection of Ha emission in young objects is com- 
monly regarded as indicative of an ongoing accretion pro- 
cess, and/or a measure of chromospheric activity. As ex- 
plained in Sect. 13.11 there is a relation between the Ha 
equivalent width EW(Hq;) and the (Ha-R) colour index 
for the objects studied bv lCNKOOl rFig.H|l. bmce our new 
very low-mass stars and brown dwarfs have similar (Ha- 
R) colours, they should also have similar Ha equivalent 
widths, with values ranging between 10 and 30 A. In TTS, 
such values are indicative of a significative accretion pro- 
cess. 

In Fig. 121 there is one remarkable outlayer: ChaHal. 
This bona-fide brown dwarf should have an equivalent 
width of ~ 20 A according to its (Ha-R) colour. However, 
it has a reported equivalent width between 34 and 99 A 
(l(^RN99l Neuhauser & Comeron [TMHI ICN KOOl -the plot- 
ted data), being by far the strongest (and most variable) 
emitter of the whole data set. Natta & Testi 1)200 l|l showed 
that Cha Ha 1 is surrounded by a disk with similar char- 
acteristics to the ones around the more massive TTS. 
However, they also report a disk around Cha Ha 2 and 
Cha Ha 9, whose Ha emission does not show a strong de- 
viation from the general correlation in Fig. (31 On the other 
hand, Cha Ha 1 also has strong X-ray emission (Neuhauser 
& Comeron [TSnSI I probably an indicator of coronal mag- 
netic activity. Therefore, its "abnormal" Ha emission 
could also be related to magnetic (chromospheric) activ- 
ity, or be the result of the combination of both accretion 
and chromospheric activity. A distinction between the two 
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processes possibly leading to the strong Ha emission, how- 
ever, cannot be made with the available data. 

Recent observations show that some objects with pho- 
tometric variability are characterized by a large Ha equiv- 
alent width (Eisloffel & Scholz^^Q. However, no photo- 
metric variability of ChaHal has been reported so far, 
although observations from different epochs are available 
(Neu hauser & Comeron nM?! I( JKJN 991 1( :JN KOOl Neuhiiuser 
et al. I2()()2|l . A study of the variability of some previously 
known brown dwarfs in Chamaeleon I was recently per- 
formed by Joergens et al. ( 20031), but these authors did not 
consider ChaHa 1 due to its poor S/N in their images. On 
the other hand, they find quite long photometric periods 
for ChaIHa2 and ChalHaS (3.2 and 2.2 days, respec- 
tively), while no variability was detected for ChaIHa4, 5, 
8 and 12. All these objects are included in Fig. 13 Since 
Cha Ha 1 is the only one showing a different behaviour in 
its Ha emission, it seems unlikely that it is connected to 
(long) photometric variability. 

Fig. 1131 shows the measured (Ha-R) colours for all our 
objects with estimated spectral type. Although no clear 
systematic behaviour is observed, interestingly, many of 
the objects with (Ha-R) colours smaller than -1 (what 
roughly corresponds to equivalent widths greater than 
30 A) have spectral types later than M7. Thus, they are 
probably brown dwarfs. Despite the larger photometric er- 
rors towards later spectral types, i.e. fainter objects, it is 
remarkable that almost all the objects with spectral type 
later than about M8 are found to have such extreme (Ha- 
R) colours. However, the low number of objects does not 
allow us to state firmly whether this is a typical property 
of substellar objects. Given that they are just a small sub- 
set of our sample (14%), it may be that these objects were 
observed near a maximum of activity. It should be noted 
that this result may simply be a consequence of the sensi- 
tivity limit of our survey. Low luminosity objects showing 
moderate to low Ha emission may appear below the com- 
pleteness magnitude and close to the detection limit of 
the Ha images, while on the contrary objects with strong 
emission are detected more easily. 

Very strong Ha emission has also been reported for 
a brown dwarf candidate in Corona Australis (Fernandez 
& Comer6n .2001.) . This object, LS-RCrA 1 (M6.5), dis- 
plays a very particular spectrum with spectral features 
of accretion and mass loss. This could also be the case 
of our strong Ha emitters, as they probably have disks 
(see Sect. 16.21 below). Such a result would be surpris- 
ing because, if brown dwarfs form from the collapse of 
a molecular cloud core, they are expected to have very 
low accretion rates. Otherwise, they would eventually ac- 
crete enough mass to start the hydrogen burning process. 
This is in contradiction with the observation of strong Ha 
emission. We note, however, that this emission can seem 
stronger in late-type objects simply due to a lower con- 
tinuum. Moreover, indications of accretion have not been 
found in all substellar objects showing a prominent Ha 
emission line. This is the case e.g. for S Ori 71, a substel- 
lar candidate member of the a Orionis cluster (Barrado y 
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Fig. 14. Upper panel: Measured Ha pseudoflux versus spec- 
tral type for our objects in Chamaeleon I. Symbols as in 
Fig. |5| The flux clearly decreases with later spectral type. 
Lower panel: logLna/Ltoi versus spectral type for our objects 
in Chamaeleon I. 



Navascues et al. 12002(1 . These authors discuss other possi- 
ble causes of the Ha emission, such as chromospheric ac- 
tivity or mass exchange between the components of a bi- 
nary system. Complementary spectroscopic observations 
of our objects will certainly help to clarify the origin and 
characteristics of their Ha emission. 

We also analysed the behaviour of the Ha flux with the 
spectral type. Since no standards were available for the Ha 
filter, we computed a "pseudoflux" using the relation: 



Fho - FHa/Fo = IQ-^H^/^.s^ 



(8) 



where mfja is the measured magnitude in the Ha filter. 
The difference between the logarithm of the pseudoflux 
Fhq and that of a hypothetical absolute flux Fhq is then 
a constant value, log Fq, depending on the reference pho- 
tometric system. 

The upper panel in Fig. E| shows log Fho versus M 
subspectral type for Chamaeleon I. A roughly linear de- 
crease of the flux with later spectral type is clearly seen, 
consistent with the progressive faintness of the objects. 
There is no indication for a change at the stellar/substellar 
boundary. The apparent saturation for the latest spectral 
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Ha-R 

Fig. 15. (Ha-R, ni(6.7)-m(14.3)) colour-colour diagram 
for the objects in our sample detected by ISOCAM in both 
passbands. The objects with very strong Ha emission have 
mid-infrared excess. 

types is most probably due to the sensitivity limit of our 
Ha survey. 

A decrease of the chromospheric activity with the spec- 
tral type has been reported by Gizis et al. I)2U(J0|I for older 
M- and L-type dwarfs. Also Zapatero Osorio et al. H2U02a|l 
found a similar behaviour in a sample of very low-mass 
members in the young a Orionis cluster. We note, how- 
ever, that these authors used Lnct/Lboi instead of the Ha 
flux in their study, and thus our result is not directly com- 
parable with theirs. Lna/Lboi is probably a better activity 
indicator, since it does not depend on the distance or the 
radius of the objects. 

The lower panel in Fig. is a plot of logLna/Lboi 
versus the spectral type for our young objects in 
Chamaeleon I. To obtain the bolometric luminosities, we 
have proceeded as in lUNKOOl using (R-I) intrinsic colours 
from Kenyon & Hartmann H1995fl and Zapatero Osorio et 
al. 1)1997(1 . Interestingly, we do not see a decrease in the 
Ha emission in these data. The situation reminds us of 
the results of Mokler & Stelzer (g002) , who studied the 
activity of a large sample of young very low-mass objects 
and found no decrease in Lx/Lboi down to spectral type 
M7, although a decrease in Lx was clearly seen in their 
data. 

Zapatero Osorio et al. H2Ul)2a|l noted that young ob- 
jects displayed, on average, higher Ha emission than their 
older field counterparts, the difference being higher for L- 
type objects. They explained this result as a consequence 
of the decline of magnetic activity and/or mass accretion 
with age. Such a direct comparison between the Ha lumi- 
nosities is not possible with our data, however, because of 
the lack of an absolute Ha photometric system. 

6.2. Evidence for Disks 

Fig. El is a plot of the colour index defined by the 
two ISOCAM passbands, m(6.7)-m(14.3) (computed with 



Eqs.iniandO, versus our measured (Ha-R) colour index. 
All our objects with strong Ha emission from their (Ha- 
R) colours are found to have mid-infrared excess according 
to Persi et al. (j2000,l . The detection of such an excess from 
young stellar objects is a well-known observational signa- 
ture of a circumstellar disk. Hence, Ha emission is related 
to the presence of an accretion disk around these objects. 

A similar proportion of the stars and brown dwarfs 
in our sample (around 30% and 37%, respectively) have 
a reported mid-infrared excess. On the other hand, while 
around 52% of the stars and 67% of the transition objects 
show near-infrared excess, only about 37% of the brown 
dwarfs do. Interestingly, mid- and near-infrared excess do 
not seem to be correlated: of the brown dwarfs with mid- 
infrared excess, just around 10% show also near-infrared 
excess. This percentage slightly increases for transition 
objects and low-mass stars (17% and 21%, respectively). 
However, the total proportion of objects detected with 
both kinds of excesses is only around 13%. This is in 
good agreement with the results of ICNKflfii who found 
that several of the previously known brown dwarfs and 
brown dwarf candidates in Chamaeleon I (selected from an 
Ha survey) showed mid-infrared excess, in spite of lacking 
near-infrared excesses. 

For TTS, the lack of near- infrared excess is usually 
attributed to an inner hole in the circumstellar disk. 
However, strong Ha emission is normally interpreted as 
a signature of the accretion process onto the surface of 
the object, what requires the presence of an inner disk. 
In a recent paper, Natta & Testi (2001) show that the 
observed spectral energy distribution of three of the pre- 
viously known brown dwarfs and brown dwarf candidates 
in Chamaeleon I can be explained if they are surrounded 
by accretion disks. They remark that, in the case of brown 
dwarfs, the presence of disks is difficult to infer from near- 
infrared photometry alone, where the emission is largely 
dominated by the photosphere. Also Oliveira et al. H20()2|l . 
who performed infrared observations of a sample of very 
low-mass objects in the crOrionis cluster, conclude that 
only a very small proportion of them show significant 
near-infrared excess. Despite this negative result, some 
of the low-mass objects in this cluster show Ha emis- 
sion and other spectral features characteristic of an ac- 
cretion process (Zapatero-Osorio et al. I2002b)l . Indeed, a 
recent L-band study has increased the disk frequency of 
the low-mass crOrionis members up to about 50%, remark- 
ably larger than inferred from the K-hand excess alone 
(Oliveira et al. lM)n|l . 

In any case, the indications for the presence of disks 
around many of our brown dwarfs and brown dwarf can- 
didates favours the hypothesis that they formed in a sim- 
ilar way as TTS, and leave open the possibility of brown 
dwarfs forming planetary systems. 

Although mid-infrared photometry is available for two 
of our faintest objects (Chal 607 and Chal 441), neither 
of them shows an excess according to the ISOCAM ob- 
servations. However, Chal 607, and also Chal 428, have 
a near-infrared excess. As argued in Sects. rOl^ and rOI 
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Fig. 16. Histogram showing the distribution of spectral 
types for the identified low-mass stars and brown dwarfs 
in Chamaeleon I. The dotted line is the distribution for 
the objects studied by previous authors (Gauvin & Strom 
Lawson et al. lTTMIirTnN99l and [nNE0(l . The local 
minimum at M2-M3 may be due to the incompleteness of 
our sample in this range (see text). 

both Chal 607 and Chal 428 may be very reddened stars. 
Their observed excess would then be due to their high ex- 
tinction (with values in the J band of 4.9 and 3.9, respec- 
tively, according to Gomez & Mardones I2U(J3|I . It seems 
thus probable that also Chal 441 is an extincted star, 
rather than a low-mass brown dwarf. 

7. The Mass Function in Chamaeleon I 

The task of constructing the substellar initial mass func- 
tion is not an easy one because the theoretical models 
still present large uncertainties for such very young ob- 
jects (< 10 Myr). Especially relevant is the choice of ini- 
tial conditions: The picture of non-accreting objects con- 
tracting from large initial radii, as used in most of the 
models, is an idealized description, which can drastically 
affect the fundamental properties (luminosity and effec- 
tive temperature) of objects younger than about 10 Myr 
(Baraffe et al. l2UU2f) . From the observational point of view, 
the main problems are the uncertainties in the luminosi- 
ties and spectral types for these young, often highly ex- 
tincted objects. Moreover, recent observations show that 
very young brown dwarfs can be variable (e.g. Eisloffel & 
Scholz l?()()2|l . Therefore, any attempt to infer an age or a 
mass from observable quantities for these objects must be 
considered with caution. 

On the other hand, at these very early stages of evolu- 
tion there is a good correlation between the spectral type 
(i.e. the effective temperature) and the mass. Thus, as a 
first step, we can study the spectral type distribution of 
the Chamaeleon I members shown in Fig. ^1 The his- 
togram includes all the new members with a determined 
spectral type from Table0] (except the few dubious objects 
discussed in the previous sections), as well as the objects 
in the surveyed area with identified spectral type from 



Ga uvin fc S trom l(TM^ . Lawson et al. l(TM5)l . ICRN99I 
and lCNKOOl The dotted line indicates the histogram for 
these previously known stars. Only the objects with spec- 
tral types between KO and L2 are shown. A numeric code 
is used: The numbers 2 to 9 correspond to KO to K7 ob- 
jects, respectively; 10-19 indicate spectral types from MO 
to M9; and numbers later than 20 are asigned to the L- 
type objects. Few objects of spectral types earlier than K 
are known in the cloud: Gauvin & Strom report four ob- 
jects of spectral type G, as well as the intermediate-mass 
stars HD 97048 (B9) and HD 97300 (AO). 

Two peaks are seen in the spectral type distribution. 
The first one corresponds to late K objects. Nonetheles, 
the decrease in the number of objects for the spectral types 
M0-M3 might not be a real characteristic of the region, 
but a consequence of the incompleteness of our sample 
in this range. We recall that the (M855-M915) colour in- 
dex saturates for spectral types earlier than M4; hence, 
twelve objects in our survey that have a spectral type 
earlier than M4 but lack a previous spectroscopic classi- 
fication are not included in Fig. E| The second peak is 
seen for spectral types around M5, approximately at the 
star/brown dwarf transition, indicating a decrease of the 
number of objects in the substellar regime. Taking M5 as 
the stellar/substellar boundary, we count around 100 stars 
(including the ones from our study whose spectral type is 
unknown) and about 51 brown dwarfs and brown dwarf 
candidates. Hence, the number of stars in the surveyed 
area is roughly doubling that of brown dwarfs. It must be 
noted, however, that this is just a lower limit to the ac- 
tual number of brown dwarfs in the cloud, because some 
low-mass objects without (detectable) Ha emission will be 
missed due to our selection criteria. Moreover, extincted 
brown dwarfs are more difficult to detect than stars, es- 
pecially at optical wavelengths. 

Taking into account the uncertainties in our abso- 
lute photometry and spectral types (recall discussions in 
Sect. l2.3l and l3.2.5|l . especially for objects earlier than M4, 
we can now construct the mass function in Chamaeleon 
and give an estimation for the value of the exponent a of 
the usual approximation for the mass function: 



dN/dM - M-" 



(9) 



For this calculation, we considered all known members 
between spectral types K7 and M7, for which our survey 
should be complete. According to the models of Baraffe et 
al. fl998), this range should correspond to masses of 1.2 
to 0.025 Mq at the age of 1 Myr. We splitted the objects 
into two bins of stars (spectral types K7-M4, with masses 
of about 1.2-0. 075 Mq) and substellar objects (spectral 
types M5-M7, with masses of about 0.075-0. 025 Mq). We 
did not take smaller bins because of the uncertainty in the 
classification for some of our objects with spectral types 
earlier than M4. By selecting two bins, we can include all 
the identified Chamaeleon I members in the calculation, 
even though we do not know their spectral types very 
accurately. There are 88 objects in the first and 39 in the 
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second bin.® With these data, we derive an exponent of 
the mass function of 

a = 0.6 ±0.1 (10) 

As mentioned above, some cloud members may have 
been missed because they do not have (measurable) Ha 
emission. To account for this possibility, we repeated the 
calculation several times including all the candidates from 
our R and I photometry that would belong to the first, 
second and both bins if they were indeed cloud members. 
These considerations lead us to the quoted error of 0.1 for 
the index of the mass function. 

Our obtained index agrees surprisingly well with the 
value obtained by Moraux et al. (2003) for the Pleiades 
cluster (0.60 ±0.11). As these authors point out, a similar 
value has been derived also for other clusters of similar 
age, thus older than Chamaeleon I. These results are also 
consistent with the values of 0.8 ± 0.4 and 0.7 ± 0.2 which 
Bejar et al. pM|l and Tej et al. found for the 

similarly old clusters (rOrionis and IC348, respectively. 
Hence, the shape of the IMF at the substellar boundary 
seems to be similar in all these young regions, regardless 
of their age and crowdedness. 

A previous estimate of the index of the substellar IMF 
in Chamaeleon I had been done bv lCNKO(H They obtained 
a value a = 1.1, larger than ours, implying a steeper in- 
crease of the mass function towards lower masses. The 
discrepancy might be due to the low number of objects 
in their study. We recall that the area covered bv KJNKdOl 
is just 15% of our whole survey, and that they studied a 
region particularly dense in very low-mass objects. 

8. Constraints on Formation Models 

The spatial distribution of the brown dwarfs in 
Chamaeleon I is found to be a natural continuation of 
the T Tauri stars towards lower masses (see Sect. 
Stellar and substellar objects also appear to have similar 
Ha emission properties fSect. lHj|l . These results, together 
with the presence of accretion disks around some brown 
dwarfs (Sect. 16.21) and the low number of possible star- 
brown dwarf systems (Sect. 0), lead us to assume that 
they most probably formed in the same way as TTS (that 
is, from the gravitational collapse of a molecular cloud) 
rather than in a circumstellar disk like planets do. 

Moreover, the fact that the low-mass objects tend 
to be concentrated in particular areas seems to indi- 
cate that they formed very near the places where we 
currently see them. We note that a star with a tan- 
gential velocity of 1 kms~^ would have moved 3 pc in 
3 Myr, which corresponds to about 1?1 at the distance of 
Chamaeleon I. Thus, if most of the brown dwarfs had been 
stellar embryos ejected from the stellar systems where they 
formed with an initial velocity of 1 kms^^, at the age of 
Chamaeleon I we should see them randomly distributed 

* We recall that the twelve objects with uncertain spectral 
type are not plotted in Fig. 1161 



over the whole surveyed region, which is not the case. The 
typical velocity of our observed young stars and brown 
dwarfs cannot be higher than 0.2-0.4 kms~^ at most. This 
value is also remarkably lower than the average ejection 
velocities of 2 kms"'^ predicted by the simulations of Bate 
et al. H2003|l . Joergens & Guenther H2001|l report a ra- 
dial velocity dispersion of 2 kms~^ for several previously 
known brown dwarf and brown dwarf candidates around 
HD 97048. They argue, however, that this value is sig- 
nificantly smaller than the radial velocity dispersion of 
the TTS in Chamaeleon I (3.6 kms~^). On the basis of 
these considerations, it seems very unlikely that the brown 
dwarfs in Chamaeleon I have been ejected from multiple 
systems. 

From the observed spatial distribution one may think 
that the existence of low-mass objects in Chamaeleon I 
could be somehow connected with the presence of the two 
intermediate mass stars. Dissipation by stellar winds from 
a nearby hot massive star could prevent the further growth 
of a stellar core, which would then remain a failed star -a 
brown dwarf. This mechanism has been invoked recently 
to explain the high concentration of substellar objects in 
Orion (Lucas & Roche IMIUl Lucas et al. lM7T|l . The stars 
in Chamaeleon I, however, are not so hot and massive as 
the ones in Orion. 

HD 97048, the one exhibiting a higher concentration 
of very low-mass neighbours, is a Herbig Ae/Be star, an 
object characterized by its strong winds. But it is doubt- 
ful whether winds from Herbig Ae/Be stars are power- 
ful enough to evaporate the envelopes of surrounding star 
forming cores. This seems even more unlikely to be cur- 
rently happening in the case of HD 97300, which might 
have reached the zero-age main sequence (e.g. Jones et al. 
I1985|) . Indeed, according to Gauvin & Strom H1992|l . ob- 
servations suggest that this part of the cloud is in a more 
evolved stage than the one containing the other interme- 
diate mass star. This is also supported by the fact that the 
cloud density, as inferred from CO measurements, is sig- 
nificantly lower in this core than in the other two (Mizuno 
et al. 1999,) . 

Furthermore, another remarkable concentration of low- 
mass objects is seen at the boundaries of the third cloud 
core, which does not contain any intermediate mass star. 
This core has not been as well studied as the other two. 
It is possible that we just observe the objects located in 
the outer part of this core, and that objects deep inside 
it are too obscured to be detected by our survey. Deeper 
observations are needed to clarify this point. If embedded 
brown dwarfs were detected, it would certainly mean that 
environmental conditions other than the presence of more 
massive stars are responsible for the formation of the sub- 
stellar objects in Chamaeleon I. 

Finally, our results from Sect. [7| indicate that the mass 
function of Chamaeleon I is rising towards lower masses. 
Hence, no matter how the brown dwarfs actually formed, 
this mechanism seems to be very efficient in producing 
substellar objects. The similar index a obtained for the 
mass functions of several young clusters and star forming 
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regions hints to a similar formation process in all these 
regions. 

9. Conclusions 

We performed a deep WFI survey of about 1.2 D° of the 
Chamaeleon I cloud to identify new very low-mass mem- 
bers of this star-forming region. Candidates were selected 
from R and I photometry and their youth was tested by 
means of their Ha-R colours. We showed that this colour 
is a useful estimator of the Ha emission rate. 

We developed a method to classify mid- to late M and 
early L-type objects from photometric observations in the 
two intermediate-band filters M855 and M915. The colour 
defined by these two passbands, (M855-M915), is corre- 
lated with the depth of the TiO and VO spectral features 
in these objects. In this way, we could provide the first 
spectral classification for about 80% of the objects in our 
sample. 

From our results, we conclude that the Chamaeleon I 
dark cloud contains a substantial amount of very low-mass 
objects, with its spectral type distribution showing a peak 
at about the star/brown dwarf transition boundary. Both 
kinds of objects have a similar spatial distribution, being 
placed predominantly near the cloud cores. Such a dis- 
tribution does not seem to be a consequence of the pres- 
ence of the two intermediate-mass stars, HD 97048 and 
HD 97300, because a large number of cloud members are 
probably associated with a dense core devoid of optical 
objects. We derive an index a = 0.6 ± 0.1 of the mass 
function in the range from 1.2 to 0.025 Mq, very similar 
to the values obtained by other authors in other young 
clusters and star forming regions. 

Both the observed spatial distribution and the scarce 
number of visual pairs among our objects seem to indicate 
that the brown dwarfs in Chamaeleon 1 formed in a way 
similar to low-mass stars, and that they were not ejected 
from their parental systems, at least with high escape ve- 
locities (greater than 0.2-0.4 kms~^). Moreover, the de- 
tection of Ha emission connected to mid-infrared excess 
shows that at least some brown dwarfs probably have ac- 
cretion disks, providing a further argument in favour of 
a star-like formation. On the other hand, near-infrared 
excess does not appear as a good diagnostic tool to the 
presence of disks around very low-mass objects. 

The Ha emitting properties of the stars and brown 
dwarfs are very similar. There is no evidence for a change 
across the substellar transition. Contrary to other authors, 
we do not see a decrease in logLna/Lboi for spectral types 
M4-M9, meaning that the latest-type objects in this range 
could have at least as strong Ha emission as earlier-type 
objects. All this hints to a common origin for the emis- 
sion in both stars and brown dwarfs, probably accretion 
processes, since both kinds of objects have disks. If it is 
confirmed, this result would contradict the commonly ac- 
cepted picture that objects with lower mass should have 
lower accretion rates (hence, weaker Ha emission) than 
those with higher masses. However, with the available data 



we cannot exclude that other mechanisms, such as chro- 
mospheric activity, are (at least partially) responsible for 
this emission. In this context, it is interesting that a sim- 
ilar result has been obtained by Mokler & Stelzer H2002|l 
for the X-ray emission of substellar objects. 
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